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A  low  heat  loss  thermal  energy  storage  system  capable 
of  storage  in  a  range  of  250-1000°F  is  described  in  this 
dissertation.   The  design  criteria  called  for  a  storage 
system  that  would  be  simple  to  construct,  use  inexpensive 
and  readily  available  materials,  utilize  an  insulating 
material  with  a  thermal  conductivity  between  0.003  -  0.015 
Btu/hrft°F  for  temperatures  ranging  between  250-1000°F,  and 
be  suitable  for  residential  or  commercial  applications. 

A  design  selected  to  fulfill  the  above  criteria  consists 
of  a  double  tank  wherein  the  annular  space  between  the  tanks 
is  filled  with  a  porous  powder  (silica)  and  then  evacuated 
to  a  pressure  of  approximately  1  mm  Hg .   Thermal  conductivity 
values  of  0.0029  -  0.110  Btu/hrft°F  were  achieved  experi- 
mentally for  vacuum  levels  ranging  from  0.5  -  3.0  mm  Hg  and 
hot  surface  temperatures  ranging  from  340  -  711°F. 


The  results  of  an  analysis  of  the  modes  of  heat  trans- 
fer through  the  porous  powder  under  vacuum  indicate  several 
means  of  potential  improvement  in  the  conductivity  of  an 
insulating  powder.   Based  on  theoretical  considerations, 
the  means  for  reducing  the  conductivity  include 


1.  selecting  powder  materials  with  a  low  thermal 
conductivity,  such  as  quartz  particles, 

2.  selecting  a  particle  size  where  an  average  void 
dimension  is  obtained  that  is  less  than  or  equal 
to  the  mean  free  path  of  the  gas  molecules  at  a 
selected  operating  condition, 

3.  selecting  gases  whose  properties  inhibit  free- 
molecular  energy  transport,  such  as  gases  with 
high  molecular  weights,  high  specific  heat 
ratios  and  large  molecular  cross-sections; 
several  of  the  more  practical  gases  being  air, 
argon  and  carbon  dioxide. 


The  results  of  this  study  indicate  for  the  temperature 
range  investigated,  250-1000 °F,  thermal  energy  storage  is 
technically  and  economically  feasible.   The  availability 
of  inexpensive  thermal  energy  storage  will  make  it  possible 
to  use  current  energy  resources  more  effectively,  particu- 
larly those  that  occur  on  an  intermittent  basis,  such  as 
direct  solar  (thermal) .   The  concept  of  thermal  energy 
storage  in  the  temperature  range  indicated  suggests  the 
possibility  of  an  independent  energy  system,  particularly 
for  residential  and  commercial  applications. 


CHAPTER  I 
INTRODUCTION 


The  more  successful  techniques  used  to  store  energy, 
electrical  or  thermal,  involve  chemical  or  mechanical  tech- 
niques such  as  lead-acid  batteries  and  thermal  energy  stor- 
age, respectively.   There  has  been  little  practical  advance- 
ment in  energy  storage  technology  over  the  techniques 
just  mentioned.   With  the  current  emphasis  on  energy 
efficiency,  there  is  a  need  to  examine  the  benefits  of 
developing  and  using  energy  storage  devices.   In  many  appli- 
cations the  need  for  thermal  energy  is  greater  than  the 
need  for  electrical  energy.   For  example,  in  a  household, 
about  85%  of  the  energy  needs  can  be  provided  by  thermal 
energy.   This  quantity  of  thermal  energy  is  needed  for 
house  heating,  house  cooling  by  absorption  air-conditioning, 
domestic  water  heating,  and  cooking.   The  remaining  15%  of 
the  energy  requirements,  for  current  lifestyles,  is  needed 
in  the  form  of  electricity  for  lighting,  refrigeration 
and  powering  communications  systems,  fans  and  pumps. 


It  should  be  noted  that  a  gas-fired  or  hot  water  driven 
absorption  air  conditioner  generally  requires  electricity 
to  run  pumps  and  controls. 


It  should  be  recognized  in  the  above  example  that 
energy,  in  the  form  of  either  heat  or  work,  is  most  effect- 
ive when  used  in  a  direct  manner.   For  example,  when  thermal 
energy  is  converted  into  electrical  energy  (work) ,  there 
follows  a  substantial  loss  as  a  result  of  the  second  law 
of  thermodynamics  in  the  form  of  the  Carnot  principle.   And 
while  electric  energy  (work)  can  be  converted  into  thermal 
energy  with  no  loss,  the  Carnot  loss  suffered  in  creating 
work  from  heat  can  be  substantial.   However,  from  a  systems 
point  of  view,  one  may  wish  to  convert  from  one  form  of 
energy  to  another  for  convenience,  to  enhance  the  system 
usefulness,  or  to  enhance  the  life-cycle  cost  of  the  system. 

Research  and  development  of  materials  for  storing 
thermal  energy  indicate  that  energy  storage  on  a  residential, 
commercial  and  industrial  scale  is  feasible.   Phase  change 
materials,  such  as  nitrate  salts,  can  store  large  quantities 
of  latent  thermal  energy  associated  with  its  phase  change  and 
organic  oils,  such  as  cottonseed  oil,  can  be  used  to  store 
large  quantities  of  sensible  thermal  energy. 

A  primary  difficulty  of  storing  thermal  energy  at 
temperatures  of  250-1000°F  is  that  one  must  find  an  econo- 
mical means  of  constructing  a  storage  system  with  low  heat 
loss  characteristics.   Very  little  literature  exists  in 
this  area.   Most  of  the  research  efforts  in  the  past 
have  been  limited  to  developing  insulation  and  insulating 
techniques  in  the  cryogenic  field. 


Many  of  the  problems  encountered  in  insulating  a 
cryogenic  vessel  are  also  encountered  in  insulating  a 
high  temperature  system.   Therefore,  the  cryogenic  area 
seems  to  be  a  good  starting  point  for  developing  a  high 
temperature  thermal  energy  storage  system.   Shared  problems 
in  the  development  of  high  temperature  or  cryogenic  storage 
container  are 


1)  heat  leaks  due  to  supports  and  tank 
penetrations , 

2)  insulation  and  tank  materials, 

3)  outgassing  of  materials  in  the  evacuated 
space, 

4)  inner  tank  support  structure, 

5)  leak-free  penetration  ports  for 
piping  and  electrical  wiring, 

6)  improved  evacuating  techniques  when 
fine  powders  are  used  as  the  insulating 
material,  and 

7)  reasonable  material  and  system  cost. 


The  solutions  to  the  above  problems  are  different  for 
cryogenic  than  for  high  temperature.   For  example,  in 
cryogenic  vessels,  some  plastics  can  be  utilized  as  a 
spacing  material;  whereas,  this  same  material  in  moderate 
or  high  temperature  environments  would  cause  problems  like 
loss  of  structural  integrity  or  outgassing. 

In  constructing  a  high  temperature  storage  vessel 
there  are  different  concepts  that  one  may  employ  to  achieve 
the  same  results.   For  example,  one  could  build  a  storage 
vessel  using  a  fibrous  insulating  material  under  a  vacuum 


but,  obviously,  this  material  would  not  support  the  inner 
vessel.   A  porous  powder  under  vacuum  could  be  used  to 
support  the  inner  vessel  and  to  insulate  the  inner  tank 
as  well.   Both  insulating  materials  may  result  in  a  similar 
heat  loss  characteristic;  however,  there  may  be  advantages 
and  disadvantages  with  using  one  material  in  place  of  the 
other. 

This  study  investigated  the  use  of  powder  insulation 
at  a  typical  pressure  level  of  1  mm  Hg  for  insulating  a 
high  temperature  thermal  energy  storage  system.   This 
technique  of  insulating  was  chosen  over  multi-layer 
insulation  because 


1)  powder  insulation  works  well  at  pressures  of 
1-10  mm  Hg  and  this  level  is  easier  to  maintain 
than  a  pressure  level  of  1  urn  Hg  or  less,  as 
required  by  multi-layer  insulation, 

2)  suitable  powders  for  withstanding  high  temperatures 
are  readily  available, 

3)  a  need  exists  for  many  uses  in  high  temperature 
technology, 

4)  it  may  be  possible  to  form  a  load  bearing  member 
of  powder  to  support  the  inner  vessel  and  this 
would  avoid  the  structural  complexities  of  build- 
ing thick  walled  vessels. 


Powder  insulations  suitable  for  high  temperatures  are 
readily  available  and  can  be  anything  from  wood  ash  to 
sand.   The  particular  material  used  in  this  study  was 
diatomaceous  earth,  a  material  readily  available,  inexpen- 
sive, and  suitable  for  high  temperature  use  because  of  low 


particle  thermal  conductivity  and  a  softening  point   of 
1600°F. 


Projected  applications  of  high  temperature  storage  systems 
The  intention  of  developing  a  high  temperature  thermal 
energy  storage  system  was  to  provide  this  component  in  a 
self-contained  cooking  system  for  residential  use  [1].   This 
system  utilizes  a  concentrating  collector  with  an  automatic 
solar  powered  tracking  mechanism  [2],  a  heat  storage  vessel, 
and  a  cooking  range.   The  concentrating  collector  coolant 
and  the  heat  storage  media  are  organic  oil. 

Four  organic  oils  (cottonseed,  peanut,  corn  and  soy- 
bean oil)  were  evaluated  for  the  above  use  [3].   In  opera- 
tion, the  organic  oil  was  heated  as  it  circulated  through 
a  heat  exchanger  located  in  the  focal  point  of  the  collector. 
With  a  minimum  temperature  of  400°F,  the  heated  fluid 
contained  in  the  storage  vessel  would  circulate  through 
the  heating  coils  of  the  cooking  range.   The  main  advantage 
of  this  system  is  that  thermal  energy  is  accumulated  and 
stored  for  use  at  a  later  time,  eliminating  the  disadvantage 
of  having  to  wait  for  direct  radiation  and  the  necessity 
to  cook  outdoors. 

The  system  can  also  be  used  as  the  basis  for  an  integra- 
ted, residential  power  system  as  shown  in  Figure  1.   For 


A  softening  point  represents  a  temperature  beyond  which 
the  yield  stress  of  a  material  drops  off  rapidly  with 
increasing  temperature. 
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example,  energy  accumulated  in  the  storage  vessel  can  be 
used  to  drive  a  heat  engine  to  produce  electricity.   The 
energy  rejected  from  the  continuously  operating  heat  engine 
can  be  used  to  drive  either  an  absorption  air  conditioner 
or  heater,  domestic  water  heater,  or  similar  low  tempera- 
ture devices.   The  high  temperature  storage  medium,  such 
as  organic  oil,  can  be  used  to  power  a  cooking  range. 

Energy  for  the  above  system  may  be  collected  from 
available  energy  resources.   Some  of  these  possible  re- 
sources are 

1)  direct  solar,  thermal  conversion, 

2)  fossil  or  nuclear  fuel  in  the  form  of  off- 
peak  electricity, 

3)  direct  solar,  electrical  to  thermal  energy 
conversion  (photovoltaic  cells) , 

4)  wind  energy  (mechanical  to  thermal  energy 
conversion) . 

As  emphasized  earlier  energy  is  best  used  in  a  direct 
manner;  that  is,  work  is  best  used  as  work  and  thermal  energy 
as  thermal  energy.   Processes  that  convert  thermal  energy  to 
work,  as  in  the  conversion  of  solar  energy  to  electric  in  a 
photovoltaic  cell  (3)  or  to  wind  (4) ,  processes  that  trans- 
form work  into  thermal  energy,  such  as  in  statements  3  and 
4  above  will  have  irreversible  losses  associated  with  them 
and  these  losses  will  be  substantial.   This  concept  is 
embodied  in  the  second  law  of  thermodynamics  and  may  be  re- 
viewed in  elementary  texts  on  thermodynamics  [4].   However,  it 
could  nevertheless  be  economically  practical  to  convert  work 


into  thermal  energy.   For  example,  if  an  industrial  process 
requires  high  temperature  fluids,  it  might  prove  advantageous 
to  purchase  electrical  energy  at  reduced  rates  from  a  utility 
company  during  periods  of  low  demand  and  to  convert  that 
power  into  thermal  energy  to  be  stored  for  later  use. 

As  apparent  from  the  above  paragraphs,  the  listed 
energy  resources  occur  on  an  intermittent  basis  and  may 
not  coincide  with  the  energy  use  pattern  for  which  they 
were  intended.   Storing  energy  for  later  use  transforms  an 
intermittent  source  of  energy  into  a  potentially  continuous 
commodity. 

In  addition  to  residential  applications  an  efficient 
cost  effective  insulating  technique  is  needed  for  high 
temperature  batteries  [5].   For  example,  a  lithium-sulfur 
battery  operates  between  707  -  797°F  and  a  sodium-sulfur 
battery  operates  between  446  -  572 °F.   A  particular  problem 
associated  with  operating  these  batteries  is  that  energy 
must  be  supplied  to  bring  them  to  operating  temperature 
and  to  maintain  stand-by  losses.   This  energy  is  not  avail- 
able for  work  and  must  be  viewed  as  lost  work.   If  an 
efficient  and  cost  effective  insulating  jacket  reduced 
these  losses,  then  high  temperature  batteries  would  be  one 
step  closer  to  reality. 

Scope  of  investigation.   The  goal  of  this  research  was  to 
develop  a  thermal  energy  storage  system  suitable  for  storing 
latent  or  sensible  thermal  energy  between  temperatures  of 
250  to  1000°F.   The  particular  objectives  required  to  ob- 
tain this  goal  were  to 


1)  find  a  practical  insulating  technique  capable 
of  reducing  the  thermal  conductivity  to  less 
than  0.015  Btu/hrft°F  @  1000°F, 

2)  identify  and  experimentally  determine  the 
necessary  thermo-physical  properties  of  the 
insulating  material, 

3)  construct  a  high  temperature  thermal  energy 
storage  system  to  verify  the  thermodynamic 
analysis  and  results  of  the  insulating  material 
and  to  verify  the  system  design. 


CHAPTER  II 
LITERATURE  REVIEW 

One  of  the  more  important  facts  to  emerge  from  the 
literature  review  was  that  no  significant  references  were 
found  concerning  efficient,  high  temperature  insulation 
suitable  for  residential  or  commercial  sized  thermal 
energy  storage  systems.   This  search  included  a  manual 
search  of  the  library  at  the  University  of  Florida 
and  a  computer  search  that  included  the  National  Aeronautics 
and  Space  Administration  (NASA)  and  the  National  Technical 
Information  Service  (NTIS)  engineering  files. 

In  view  of  the  above,  the  goal  of  this  research  was 
to  develop  a  thermal  energy  storage  system  suitable  for 
storing  latent  or  sensible  thermal  energy  between  tempera- 
tures 250  -  1000°F.   The  successful  development  of  a 
thermal  energy  storage  system  depends  upon  the  type  and 
quality  of  the  insulation  used  in  constructing  the  system. 
The  initial  literature  review  shows  that  there  are  basically 
three  types  of  insulation  suitable  for  the  indicated 
temperature  range.   Based  on  physical  structure  they  can 
be  described  as 

1)  fibrous  insulation, 

2)  laminated  insulation, 

3)  granular  insulation. 
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Fibrous  insulation  includes  such  materials  as  fiberboard 
and  rockwool,  which  are  suitable  for  low  and  moderate 
temperature  applications.   Laminated  insulation  consists 
of  alternating  layers  of  low  density  insulation  such  as 
ceramic  paper  and  low  emittance  radiation  shields.   Granu- 
lar insulation  includes  such  materials  as  sand,  micro- 
sheres  of  glass,  and  diatomaceous  earth.   Laminated  and 
granular  insulation  are  suitable  for  use  from  cryogenic  to 
high  temperature  applications. 

As  described  by  Jakob  [6]  ,  some  of  the  earliest  work 
on  conduction  through  porous  material  was  done  by  such  men 
as  Maxwell,  Nusselt  and  Eucker.   These  investigators  included 
gross  simplifications  of  the  heat  transfer  mechanisms  in 
their  analyses.   Their  aim  was  to  make  it  possible  to  obtain 
reasonable  answers  to  a  very  complex  problem  with  a  minimum 
of  data.   It  seems  that  little  advancement  has  been  made  in 
the  area  of  heat  transfer  through  porous  materials  since 
these  early  studies. 

The  basic  ideas  presented  by  Reference  6  concerning 
Maxwell's,  Nusselt 's,  and  Eucker 's  analyses  indicate  that 
all  three  viewed  a  porous  body  as  either  a  base  material 
with  small  particles  of  another  substance  distributed 
throughout  that  material,  or  as  alternating  layers  of  solids 
and  gases  with  heat  flow  perpendicular  to  the  layers. 

One  of  the  classic  works  on  thermal  conductivity  of 
fine  powders  at  reduced  pressures,  as  described  by  Jakob  [6], 
was  that  of  Smoluchowski .   The  essence  of  Smoluchowski ' s 
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work  was  that  at  reduced  pressures,  gas  molecules  would 
experience  slip  flow  in  the  vicinity  of  a  solid  surface. 
This  causes  a  discernible  temperature  discontinuity  to 
exist  between  the  solid  surface  and  the  gas.   When  this 
theory  is  applied  to  fine  powders,  one  correctly  concludes 
that  many  temperature  discontinuities  are  created.   These 
added  resistnaces  in  turn  reduce  the  transfer  of  heat 
through  porous  powder. 

Most  of  the  early  works  of  heat  flow  through  porous 
material  were  developed  considering  ambient  or  lower  tempera- 
tures.  Few  studies  were  found  concerning  the  flow  of  heat 
from  one  particle  to  another.   One  study,  by  Strong  et  al. 
[7],  analyzed  heat  flow  through  glass  fibers  at  a  pressure 
of  10  pm  Hg.   Even  though  this  investigation  worked  with 
materials  with  smooth  surface  geometry,  its  experimental 
particle-to-particle  conductivity  was  off  by  an  order  of 
magnitude  when  compared  with  theoretical  values  of  heat 
flow  perpendicular  to  the  fiberglas  rods. 

In  the  first  part  of  the  twentieth  century,  two  facts 
important  in  insulation  were  discovered  and  then  investiga- 
ted.  First,  it  was  noted  by  Smoluchowski  [8]  that  soot 
had  a  lower  conductance  than  air.   At  that  time  it  was 
thought  that  the  conductivity  of  air  established  a  minimum 
value  for  thermal  insulation.   However,  Kistler  and  Caldwell 
[9]  in  1934  found  that  silica  aerogel  had  a  thermal  con- 
ductivity less  than  air  and  reasoned  that  the  lower  conduc- 
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tivity  resulted  from  the  fact  that  the  pore  size  of  the 
material  was  less  than  the  mean  free  path  of  molecules  of 
air. 

The  second  important  fact,  noted  by  Kistler  and 
described  by  Wilkes  [10],  was  that  when  commercial  cork 
was  placed  in  an  evacuated  chamber,  thermal  conduction  and 
radiant  transfer,  through  the  material,  at  low  temperature, 
was  found  to  be  a  small  part  of  the  overall  heat  transfer. 
This  indicated  that  air  was  responsible  for  the  majority  of 
the  energy  transport  through  the  material. 

In  a  more  recent  study,  Peterson  [l]  showed  that 
successive  thin  foil  radiation  shields  of  low  emissivity 
wound  with  fine  steel  wool  spacers  were  exceptionally  good 
insulators,  especially  when  used  under  vacuum  conditions. 

Based  upon  the  three  developments  mentioned  above,  it 
was  possible  to  construct  super- insulated  vessels,  at  least 
for  cryogenic  purposes.   The  requirements  for  space  age 
technology  placed  a  premium  on  lightweight,  high  efficiency 
insulation,  and,  consequently,  great  interest  was  shown  in 
cellular  and  laminated  or  multi-layer  insulation  [12,  13]. 
Most  of  the  development  work  on  multi-layer  insulation 
was  done  for  cryogenic  vessels,  but  some  work  was  performed 
for  high  temperature  applications .   For  example,  the 
early  design  of  a  thermoelectric  nuclear  power  system 
for  spacecrafts  used  multi-layer  insulation  operating 
at  a  hot  surface  temperature  of  1285°F  [14].   The  multi- 
layer consisted  of  alternating  layers  of  copper  foil  with 
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quartz  paper  separators,  and  had  an  upper  temperature  limit 
of  1400°F.   The  vacuum  space  was  designed  and  built  to 
hold  a  pressure  of  less  than  10  ym  Hg  for  at  least  five 
years . 

As  evident  from  much  of  the  literature,  multi-layer 
insulation  for  spacecraft  and  for  cryogenic  work  was 
expensive.   Prior  to  the  need  for  super-insulators  for 
spacecraft  use,  evacuated  powders  were  sufficient  for  most 
low  temperature  pruposes.   As  pointed  out  by  Black  et  al. 
[15],  although  powder  insulation  does  not  have  as  low  a 
thermal  conductivity  as  multi-J.ayer  insulation,  its  advantages 
are 

1)  it  can  operate  effectively  at  a  few  millimeters 
of  mercury  or  less,  thus  eliminating  the  need 
for  high  vacuum, 

2)  it  can  be  used  in  vessels  of  intricate  design,  and 
this  eliminates  the  wrapping  and  bonding  of  each 
component  in  a  vessel,  as  is  necessary  with  multi- 
layer insulation, 

3)  it  costs  less  than  multi-layer  insulation. 

In  view  of  this,  it  seemed  that  porous  powder  under  partial 
vacuum  should  be  suitable  for  terrestrial  high  temperature 
insulation  applications  where  weight  is  not  a  serious  dis- 
advantage . 

There  is  a  growing  need  for  efficient  high  temperature 
insulation.   The  present  technique  used  to  insulate  large, 
high  temperature  (200-1000°F)  storage  vessels  is  to  use 
thick  layers  of  granular  insulation  at  atmospheric  pressure 
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[16,17],   This  approach  seems  reasonable  for  large  thermal 
energy  storage  systems  in  view  of  the  quantity  of  energy 
stored  relative  to  the  amount  of  energy  lost.   In  smaller 
systems,  such  as  a  residential  sized  thermal  energy 
storage  vessel  capable  of  storing  106  Btu  with  minimal 
losses,  the  storage  efficiency  is  more  closely  tied  to  the 
overall  system  size  and  shape.   This  results  from  a  smaller 
volume-to-surface  ratio.   No  references  were  found  con- 
cerning high  efficiency  insulation  suited  to  small  systems. 
Available  literature  indicates  that  little  consideration 
has  been  given  to  high  temperature  thermal  energy  storage 
systems.   It  seems  that  more  attention  has  been  directed 
into  developing  heat  storage  media  than  in  developing  the 
heat  storage  media  containment  system. 

Table  1  has  been  compiled  to  illustrate  the  types  of 
materials,  and  their  associated  thermal  conductivities, 
that  are  available  as  high  temperature  insulators.   Note 
that  most  high  temperature  insulating  materials  listed  in 
Table  1  have  similar  thermal  conductivities.   This  result 
is  not  unexpected  since  these  materials  are  composed  of 
similar  materials  which  is  basically  Si02.   All  of  the 
above  materials  are  at  atmospheric  pressure,  with  the  excep- 
tion of  the  diatomaceous  earth  which  is  at  1  mm  Hg .   The 
reduced  pressure  results  in  a  thermal  conductivity  value 
that  is  approximately  one  order  of  magnitude  lower  than 
similar  material  at  atmospheric  pressure.   Lower  thermal 
conductivities  can  be  obtained  as  will  be  explained  in  the 
next  chapter. 
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CHAPTER  III 
THEORETICAL  CONSIDERATIONS  OF  POWDER  INSULATION 

Introduction 

It  is  relatively  easy  to  measure  the  overall  apparent 
conductivity  of  a  powder  insulation.   However,  if  one 
wishes  to  alter  and  optimize  apparent  conductivity  other 
than  by  increasing  thickness,  it  is  necessary  to  have 
an  acceptable  theory  with  which  to  proceed.   If  one  ideal- 
izes the  properties  of  the  media  which  form  the  powder  and 
spaces,  then  one  can  model  its  apparent  conductivity. 
Unfortunately,  the  true  geometry  is  not  well  modeled  by 
a  simple  approximation,  nor  are  the  properties  adequately 
approximated  as  ideal.   As  a  result,  it  is  very  difficult, 
perhaps  impossible  to  evolve  an  analytically  tractable 
model  that  also  accurately  describes  the  phenomena  as 
they  occur  in  the  powder.   This  seems  to  be  the  case  with 
modeling  of  energy  transport  through  porous  powder  insula- 
ting material. 

In  this  study  diatomaceous  earth  (981  Si02)  was  used 
as  the  insulating  powder.   In  view  of  the  above  discussion 
the  major  questions  encountered  in  modeling  the  irregularly 


* 

Apparent  conductivity  is  heat  flow  from  conduction,  con- 
vection, and  radiation  combined  and  used  in  equations  as 
if  it  were  conduction  alone. 
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shaped  particles  were  how  to  describe  factors  that  enhance 
or  impede  heat  flow  through  this  material.   These  factors 
are 

1)  gaseous  thermal  conduction,  as  limited  by  the 
molecular  mean  free  path  and  influenced  by  gas 
pocket  geometry  and  boundary  conditions, 

2)  conduction  through  solids,  limited  especially 
by  particle  geometry  and  boundary  conditions, 

3)  radiation  resistance  determined  by  particle- 
void  space  geometry,  particle  surface  properties 
and  boundary  conditions. 

The  best  understood  of  the  processes  applicable  to  heat 
transfer  through  powders  seems  to  be  gaseous  conduction, 
followed  by  radiation  and  particle-to-particle  conduction. 
The  relative  importance  of  each  process  is  discussed  later. 
The  remainder  of  this  chapter  is  devoted  to  developing  a 
model  for  heat  flow  through  porous  powder  insulation  which 
addresses  the  above  questions. 

Theoretical  Model 

A  mechanistic  model  for  heat  flow  through  a  porous 
powder  is  sh6wn  in  Figure  2.   The  total  energy  flow  through 
the  gas -powder  mixture  is  a  combination  of  conductive, 
convective,  and  radiant  interaction  between  the  solid  and 
gaseous  materials  present.   Symbolically,  the  total  heat 
flow  through  the  gas  powder  mixture  can  be  written  in  func- 
tional form  as 
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qT  -  qT(qc,  qr,  qh) 


(i) 


For  powders  with  particle  size  less  than  0.2  inches, 
investigators  such  as  Hill  and  Wilhelm  [20]  found  no 
evidence  of  convection  in  quiescent,  gas-solid  beds  for 
conditions  of  interest  to  this  study.   Also,  in  this  work, 
the  powder  material  was  studied  at  pressures  ranging  from 
atmospheric  to  less  than  1  mm  Hg.   Consequently,  convection 
heat  transfer  was  eliminated  from  further  consideration. 


/// 


////// 


undary 


Figure    2 

Mechanistic  Description  of  xHeat  Transfer  Through 
Porous  Powder 

(1)  Gaseous  conduction,  (2)  solid  conduction, 
(3)  radiant  heat  transfer  through  voids,  (4)  radiant 
heat  transfer,  solid- to-solid,  (5)  conduction, 
point  contact,  (6)  convection. 
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The  justification  for  eliminating  convection  heat 
transfer  under  the  conditions  discussed  above  can  be 
shown  by  estimating  the  convective  heat  transfer  coef- 
ficient, h,  between  particles.   This  analysis  is  given 
in  Appendix  B.   The  results  of  this  analysis  reveal  that 
the  Nusselt  number, 

Nu  =  *£  (2) 

is  of  order  unity  at  atmospheric  pressure  and  at  tempera- 
tures between  250  -  1000°F,  then  the  convective  heat 
transfer  coefficient  reduces  by  two  orders  of  magnitude. 
Hence,  convection  heat  transfer  at  high  temperature  and 
low  pressure  can  be  eliminated  from  further  consideration. 
In  view  of  Equation  2,  energy  transport  through  a  gas-solid 
mixture  in  which  convection  has  been  eliminated  reduces 
to  primarily  a  function  of  conduction.   It  will  be  shown 
later  that  net  radiant  transfer  contributes  only  a  small 
fraction  of  energy  transfer  under  conditions  encountered 
in  this  work.   Thus,  Equation  1  is  rewritten  as 

qT  ■■  qTCqc)  C3) 

The  remainder  of  this  chapter  is  devoted  to  examining 
in  detail  conduction  through  both  gas  and  solids  and  the 
small  contribution  of  radiant  transfer  through  porous  pow- 
ders (silica)  under  partial  vacuum. 
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Thermal  Conductivity  of  a  Perfect  Gas 

Introduction.   This  section  examines  the  thermal 
conductivity  of  a  perfect  gas.   It  is  implicit  in  the  use 
of  the  property  thermal  conductivity  that  a  continuum 
exists.   A  continuum  is  assumed  to  exist  as  long  as  the 
mean  free  path  of  a  molecule  is  comparable  to  the  smallest 
significant  dimension  of  its  surroundings.   The  ratio  of 
the  mean  free  path,  \,    of  a  gas  molecule  to  the  average 
dimension  of  its  boundaries,  d~,  is  referred  to  as  the 
Knudsen  number.   That  is, 

Kn  =  X/cT.  (4) 

This  number  is  useful  in  identifying  the  conditions  under 
which  continuum  assumptions  apply.   The  need  for  identify- 
ing whether  a  continuum  exists  is  important  in  determin- 
ing the  proper  relationship  required  for  calculating  the 
thermal  conductivity  of  a  gas.   As  will  be  shown  in  the 
following  sections,  the  thermal  conductivity  of  a  gas  in  a 
gas-powder  mixture  is  influenced  by  the  magnitude  of  the 
Knudsen  number. 

Conductivity  in  a  perfect  gas.   This  section  examines 
thermal  conductivity  of  a  continuum.   A  continuum  is  said 
to  exist  when  the  mean  free  path  of  a  gas  molecule  is  small 
(Kn  <  <  1)  compared  to  the  dimensions  of  its  boundaries. 
For  energy  transport  in  a  continuum,  the  gaseous  thermal 
conductivity  is  given  by 
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a    =  -  K  VI  (5) 

lgc       S 

The  conductivity  coefficient,  Kff,  is  given  by  [21] 

Ka  =  %(9Y  -  5)(%PCCvX)  (6) 


Density,  p,  is  given  by 


PM 

P    R  T 
u 


(7) 


where  P  is  the  pressure,  M  is  the  molecular  weight  of  the 

oas   R   is  the  universal  gas  constant,  and  T  is  the  tempera- 
&   '   u 

ture  of  the  gas  (absolute  temperature) .   The  average  gaseous 
molecular  velocity,  C,  is  given  by 

8  R  T  , 

The  mean  free  path,  X,  for  a  perfect  gas  is  given  by 
R  T 


X  = 


(9) 


/2Va2A  P 
n 


where  a  is  the  collision  diameter  of  the  gaseous  molecule, 
P  is  pressure  and  A   is  Avogadro's  number. 

ilume,  Cy,  is  given  by 

C   =  R  /M(y-1)  (10) 

v    u'  v    J 

where  Y  is  the  specific  heat  ratio. 

Combining  Equations  7,  8,  9,  and  10  into  Equation  6 
yields 


23 


g     4(Y-i)  lwJ    ^^r  (11} 

Equation  11  indicates  that  gaseous  conductivity  is  a 
function  of  temperature  and  that  it  should  rise  with  an 
increase  in  temperature  according  to  T2.   With  an  increase 
in  temperature  over  the  range  of  250  -  1000°F,  and  taking 
into  account  the  temperature  variation  of  C   and  y,    experi- 
mental values  of  thermal  conductivity  of  air  at  atmospheric 
pressure  change  from  0.0192  to  0.0337  Btu/hrft°F. 

From  Equation  11  one  can  identify  the  gaseous  molecular 
properties  that  will  cause  the  conductivity  coefficient 
to  be  small.   For  example,  one  would  choose  a  gas  with  high 
molecular  weight,  large  collision  diameter,  and  a  high 
specific  heat  ratio.   Table  2  lists  some  values  of  M,  a 
and  y  for  several  of  the  more  common  gases. 

It  is  interesting  to  note  in  Equation  11  that  the 
thermal  conductivity  of  a  perfect  gas  is  independent  of 
pressure.   Mathematically,  the  reason  for  this  result  is 
that  the  direct  proportionality  of  pressure  in  the  density 
term,  p,  cancels  with  the  inverse  proportionality  of  pres- 
sure in  the  mean  free  path,  X,    expression.   Physically, 
this  means  that  if  the  density,  or  number  of  gas  molecules 
in  a  given  space  is  decreased,  the  mean  free  path  or  dis- 
tance between  collisions  of  gas  molecules  increases  so  that 
the  transport  of  energy  is  unchanged. 
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TABLE  2 

Molecular  Weight,  Collision  Diameter,  and  Specific 
Heat  Ratio  of  Selected  Gases 


Molecule  (lb  /lb  -mole) 
m   m 

A  39.948 

Air  28.97 

CO  28.011 

C09  44.011 

H  2.015 

He  4.003 

N_  28.013 

02  31.999 

Ne  20.183 


Specific  heat  ratios  for  gases  at  low  pressure  at  80°F. 
Corrected  value  [24] . 


o    [22] 
(ft  x  10   ) 

(d 

T  [23]* 
imensionless) 

0.125 

1.667 

1.1975 

1.40 

1.05 

1.40 

1.10 

1.285 

.886 

1.40 

.689 

1.667 

1.23 

1.40 

1.94  + 

1.40 

1.61 

1.667 
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The  calculation  of  the  thermal  conductivity  of  a  gas, 
as  calculated  with  the  perfect  gas  continuum  assumptions, 
yields  results  that  are  adequate  for  describing  real 
gases  in  many  situations.   However,  once  the  variables  of 
temperature,  pressure,  and  boundary  conditions  are  such 
that  a  continuum  no  longer  exists,  that  is  Kn  >  1,  then 
Equation  6  must  be  modified.   This  is  the  topic  of  the 
following  section. 

Thermal  Conductivity  of  a  Non-Continuum  Gas 

As  mentioned  in  the  previous  section,  the  conductivity 
of  gases  is  theoretically  independent  of  pressure  in  a 
continuum.   However,  both  Knudsen  [25]  and  Smoluchowski 
[8]  developed  theories  concerning  gaseous  thermal  conduction 
for  non -continuum.   Non-continuum  and  in  the  limit  free- 
molecular  energy  transfer  occurs  when  the  confining 
boundaries  of  a  gas  molecule  become  smaller  than  the 
mean  free  path  of  the  gas  molecule.   For  this  condition 
another  phenomenon  occurs  that  limits  the  energy  transfer 
between  a  gas  molecule  and  a  solid  boundary.   Knudsen 
found  that  when  a  gas  molecule  strikes  a  solid  surface 
that  the  gas  molecule  does  not  necessarily  come  into 
thermal  equilibrium  with  the  surface.   That  is,  a  tempera- 
ture discontinuity  was  found  between  the  temperature  of  the 
surface  and  the  temperature  of  the  gas  molecule  leaving 
the  surface.   Thus,  Knudsen  introduced  a  constant  surface 
property,  called  the  thermal  accommodation  coefficient,  a, 
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which  is  a  measure  of  the  degree  of  energy  exchange  between 
a  gas  molecule  and  a  solid  surface.   The  maximum  value  of 
the  accommodation  coefficient  is  unity.   This  value  is 
approached  by  surfaces  that  are  roughened. 

The  energy  transfer  between  parallel  plates  a  distance 
x  apart  and  at  moderately  low  pressure  is  given  by  [26] 

K  (T  -T  ) 
*  =   S  +  2BX  (12) 

where  6  is 


and 

e  =  -j  (9y-5)  (14) 


At  very  low  pressures,  in  the  range  of  molecular  conduction, 
it  was  found  that  energy  transfer  becomes  independent 
of  the  plate  spacing,  and  proportional  to  the  pressure,  P. 
Rewriting  Equation  9  for  the  mean  free  path  of  a  gas 
molecule  as 

A      =   B/P  (15) 


where 


R   T 

B   = (16) 

/2rra2A 
n 
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and  rewriting  Equation  12  in  view  of  Equation  15  for  condi- 
tions where  energy  transfer  is  independent  of  x  yields, 


K  (T„-Tr)P 


Hence,  molecular  heat  conduction  is  proportional  to  pressure 
at  very  low  pressures,  (Kn  >  1) . 

Equation  17  did  not  appear  suitable  for  the  present 
investigation  for  the  following  reasons: 

1.  an  equation  was  needed  that  would  be  suitable 
for  energy  transfer  for  continuum,  transitional, 
and  non-continuum  conditions,  and  not  just  the 
latter , 

2.  the  use  of  Equation  17  requires  a  knowledge  of 
accommodation  coefficients  and  these  values  are 
generally  not  available, 

3.  Equation  17  was  derived  on  the  basis  of  parallel 
plates  and  the  material  under  investigation  is 
powder  material  with  particles  and  cavities  that 
are  irregular  in  size  and  shape  and  oriented  at 
random;  that  is,  the  boundary  dimensions  must  be 
taken  into  account. 

Strong  et  al.  [7]  have  found  that  Equation  6,  when  modi- 
fied for  non-continuum  conditions,  gave  experimental  results 
that  were  in  good  agreement  with  theoretical  predictions 


28 

over  the  range  of  continuum,  transitional,  and  non-continuum 
conditions.   The  material  investigated  by  Strong  et  al . 
[7]  was  fiberglass  rods  which  formed  small  cavities,  as  does 
the  powder  material  of  this  investigation.   The  basic 
assumptions  in  using  Equation  6  were 

1.  the  accommodation  coefficient  could  be  approxi- 
mated as  unity,  and 

2.  that  by  modifying  the  definition  of  the  mean  free 
path,  the  equation  would  be  suitable  over  the 
range  of  pressure  from  continuum  to  non-continuum 
conditions . 

Since  the  accommodation  coefficient  is  a  measure  of  the 
degree  of  energy  transfer  between  a  gas  molecule  and  a 
solid  surface,  then  it  seemed  reasonable  to  conclude 
that  for  a  material  with  cavities  on  the  order  of  the  mean 
free  path  of  the  gas,  that  an  effective  accommodation 
coefficient  near  unit  will  result.   That  is,  a  gas  mole- 
cule will  interact  with  the  boundaries  of  a  given  cavity 
a  sufficient  number  of  times  to  reach  equilibrium  before 
escaping  to  another  cavity. 

In  a  non-continuum  gas  there  are  two  space  dimensions 
that  must  be  taken  into  account  —  the  mean  free  path,  X, 
of  the  gas  molecule  and  the  average  space  or  void  dimension 
of  its  confining  boundaries,  3".   When  the  mean  free  path 
is  less  than  the  space  dimension,  d~,  Kn  <  1,  perfect  gas 
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conductivity  is  independent  of  pressure,  as  expressed 
by  Equation  6.   When  the  mean  free  path  of  the  gas  molecule 
is  close  to  or  greater  than  the  space  dimension,  Kn  >  1, 
then  the  conduction  equation  must  be  modified  to  take  this 
into  account.   A  modified  mean  free  path  of  a  gas  molecule, 
L,  is  related  to  the  continuum  mean  free  path,  A,  and 
the  gas  space  dimension,  3",  by  Strong  et  al .  [7]  as 

L  =  (r  +  =)_1  (18) 

X        d 

In  the  continuum  region  where  A  <  <  d,  then 


In  the  non-continuum  region  where  A  >  >  d",  then 

L  -  1 

To  relate  the  modified  mean  free  path,  L,  to  gaseous  con- 
ductivity, recall  that 

K   =  %(9Y-SmpCCvA)  (6) 

Substituting  Equation  12  into  Equation  6  yields 

Ka  -  %(9Y-5)  i2pCCv  (1/A  +  1/d)"1  (19) 

where  K  is  the  modified  gaseous  thermal  conductivity  coef- 
ficient.   Equation  19  can  now  be  used  to  calculate  the 
thermal  conductivity  of  a  gas  for  continuum,  transitional, 
and  non-continuum  conditions. 
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Finally,  to  show  that  gaseous  thermal  conductivity, 
K  ,  is  related  to  pressure,  recall  Equations  15  and  16 
and  with  the  appropriate  substitution,  Equation  19 
becomes 


xr^TfcSfr  Cl/P**  l/B) 


K   =  [9y'SA     (-r£)h    Cl/PcT  +  1/B)"1  (20) 


Thus,  as  the  pressure  is  lowered,  at  a  fixed  void  geometry 
and  temperature,  the  gaseous  thermal  conductivity  decreases 
monotonically.   Note  that  as  the  product  of  PcT  becomes 
smaller  than  B,  the  former  value  dominates,  and  it  can  be 
concluded  that  X     varies  essentially  linearly  with  pressure. 
Conversely,  as  the  product  of  PcT  becomes  large  relative  to 
B,  K   is  seen  to  be  essentially  independent  of  pressure. 

Equation  20  also  suggests  that  the  gaseous  conductivity 
will  remain  the  same  if,  at  a  given  temperature,  the  modified 
mean  free  path,  L,  is  held  at  a  constant  value,  that  is, 

L  =  (1/P<T  +  1/B)"1  =  constant  (18) 

This  means  that  a  given  value  of  K*  can  be  maintained 
by  either  decreasing  d"  and  increasing  P,  or  by  increasing 
3"  and  decreasing  P.   This  suggests  that  by  decreasing  the 
particle  size,  which  will  lower  the  average  void  dimension, 
cF,  a  higher  gas  pressure  will  be  acceptable.   For  example, 
by  selecting  particles  such  that  the  pressure  level  is 
raised  from  0.10  mm  Hg  to  5  mm  Hg ,  while  maintaining  a 
constant  value  of  X   ,  many  of  the  problems  associated  with 
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low  pressure,  such  as  the  need  for  diffusion  pumps,  out- 
gassing  of  surfaces,  and  long  term  leakage,  can  be  avoided 
or  minimized.   It  is  the  control  over  the  particle  size, 
and  consequently  the  gas  space  dimension,  that  makes  a 
gas-powder  mixture  attractive  for  high  temperature  insula- 
tion.  Both  Equations  17  and  20  indicate  that  molecular 
energy  transport  can  be  reduced  by  selecting  gases  with 
high  molecular  weights  and  high  specific  heat  ratios. 
This  same  trend  is  also  exhibited  for  conduction  of  heat 
by  gases  in  a  continuum.   Equations  17  and  20  do  not  show 
a  dependency  on  the  collision  cross-section  in  the  regime 
of  molecular  conduction  (Pd  <  <  B)  since  the  mean  free 
path  has  been  superseded  by  the  confining  boundaries  of 
the  gas  molecule.   Also,  for  non-continuum  or  molecular 
conduction,  T     is  found  to  vary  with  temperature  by  T  2. 
Hence,  molecular  conductivity  decreases  as  the  temperature 
is  elevated. 

When  the  thermal  conductivity  of  a  gas  has  been  reduced 
below  the  particle-to-particle  thermal  conductivity  in  a 
gas  powder  bed,  little  will  be  gained  by  further  reducing 
conduction  through  the  gas.   Under  these  circumstances, 
the  particle-to-particle  thermal  conductivity  sets  the  lower 
limit  of  the  apparent  thermal  conductivity  through  a  porous 
powder  under  vacuum.   (It  will  be  shown  later  that  net 
radiant  transfer  appears  negligible  for  the  conditions 
encountered  in  this  study) . 
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The  following  section  examines  the  contribution  of 
particle-to-particle  thermal  conductivity  to  the  apparent 
thermal  conductivity  of  a  powder  under  vacuum. 

Particle-to-Particle  Conductivity  of  a  Porous  Powder 

Introduction.   Conduction  through  a  region  of  contact 
between  particles  is  not  easily  analyzed  except  for  ideal 
geometries.   For  an  idealized  geometry  one  may  express 
energy  transfer  between  particles  of  given  size,  mathema- 
tically define  the  contact  area,  and  express  the  proper 
number  of  contact  points.   In  most  real  powders,  for  example, 
see  Figure  8,  page  68,  the  particles  are  both  irregular  in 
shape  and  size  and  oriented  at  random.   These  factors  make 
it  difficult  to  determine  the  parameters  necessary  to 
express  energy  transport  by  conduction  through  the  powder. 

Even  with  the  difficulties  of  analyzing  particle-to- 
particle  conductivity  mentioned  above,  and  the  fact  that 
there  is  little  in  the  literature  analyzing  this  topic,  it 
is  useful  to  develop  an  approximate  mathematical  expression 
to  describe  particle-to-particle  conductivity.   With  a  mathe- 
mathical  expression  one  may  readily  identify  those  variables 
that  control  the  behavior  of  the  material  under  study. 

Theoretical  Model .   For  particle-to-particle  conductivity 
to  exist  there  must  be  an  unbroken  path  connecting  the  parti- 
cles and  continuing  through  them.   It  is  through  this  path 
that  energy  is  transferred  by  solid  conduction.   A  particu- 
lar problem  in  analyzing  heat  transfer  through  a  powder  that 
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consists  of  irregularly  shaped  material,  such  as  the  material 
under  consideration  in  this  study,  is  the  difficulty 
of  mathematically  defining  the  path  through  which  energy 
flows.   For  example,  in  Figures  8,  9,  and  10  (pages  68 
and  6.9)  one  sees  a  variety  of  shapes  ranging  from  per- 
forated, hollow  and  solid  cylinders  to  perforated  plates, 
and  irregularly  shaped  solid  pieces.   With  the  variety  of 
shapes  shown,  following  the  path  of  heat  transfer  from 
particle-to-particle  appears  to  be  virtually  an  impossible 
task,  thus  the  need  for  a  simplifying  model. 

In  following  the  path  of  energy  transfer  from  one 
particle  to  the  next,  one  finds  that  the  region  of  greatest 
resistance  occurs  at  the  smallest  area,  that  is,  in  a  region 
of  contact.   If  the  dimensions  of  the  contact  region  are 
small  relative  to  the  bulk  particle  size,  then  one  may 
approximate  particle-to-particle  energy  transfer  in  a 
powder  by  analyzing  the  resistance  due  to  the  contact 
region  only.   This  concept  is  best  illustrated  by  an 
example . 

Consider  a  cylinder  as  a  particle  body  as  shown  in 
Figure  3.   Beside  this  cylinder  is  a  smaller  cylinder  which 
represents  the  contact  region.   The  two  cylinders  are  in 
intimate  contact  and  they  are  aligned  on  the  same  axis. 
Each  cylinder  has  a  length  to  diameter  ratio  of  10:1  and 
the  size  of  the  two  cylinders  relative  to  each  other  is 
100:1.   Assuming  the  conductivity  coefficient  to  be  the  same 


Contact  region 


(a) 


General  Particle  Shapes  In  Contact 
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Particle  interface 


Particle  1 


Particle  2 


Contact  region 
(b) 

Idealization  Of  Particles  In  Contact 

Figure  3 

Idealization  Of  A  Contact  Region 
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for  both  cylinders,  one  can  easily  demonstrate  that  the 
resistance  to  heat  flow  in  the  smaller  cylinder  is  100 
times  that  of  the  larger  cylinder.   This  indicates  that  the 
contact  resistance,  as  represented  by  the  smaller  cylinder, 
predominates  and  that  the  body  resistance  of  the  particles, 
as  represented  by  the  larger  cylinder,  is  negligible.   In 
Chapter  V  it  will  also  be  demonstrated  that  the  relative 
dimensions  of  the  particle  body  and  contact  region  are 
such  that  a  resistance  ratio  greater  than  100:1  is  not 
unusual . 

Note  carefully  that  a  distinction  is  made  between  a 
contact  area  and  a  region  of  contact.   A  contact  area  is 
the  physical  area  or  interface  between  two  particles.   A 
contact  region  is  defined  to  be  the  contact  area  along  with 
the  particle  materials  adjacent  to  the  contact  area  that 
restrict  energy  transfer  between  two  particles  forming  the 
contact  area.   See  Figure  3.   The  contact  region  is  viewed 
as  a  continuous  path  of  solid  material  through  which  con- 
duction heat  transfer  occurs  even  though  an  interface 
exist  between  two  particles.   This  assumption  is  based  upon 
the  existence  of  sufficient  stresses  to  maintain  the  two 
particles  in  intimate  contact.   That  is, 

S  =  F/A  (21) 

where  S  is  the  average  stress  between  the  particle,  F 
is  the  force  or  weight  supported  by  a  particle,  and 
A  is  the  contact  area  between  two  particles.   It  will 
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be  demonstrated  in  Chapter  5  that  sufficient  stress  exists 
at  a  contact  region  to  assume  intimate,  solid  contact 
between  particles. 

The  contact  region  of  a  solid  particle  can  be  analyzed 
by  any  suitable  coordinate  system.   For  simplicity,  cylindri- 
cal coordinates  will  be  used.   Assuming  constant  conducti- 
vity, one  can  express  energy  conservation  in  conduction  by 

32T  +  _3T_  +   32T    32T    q '"  u    ldj_  f 

3T2"   r3r   rzez   ax"2"   Tc     <x3t  {      } 

s 

Since  the  contact  region  is  small  relative  to  the  particle 

body  dimensions,  with  the  latter  varying  from  1  ym  to 

lOOym  in  diameter,  the  particle  material  conductivity 

coefficient,  k    is  assumed  to  be  homogeneous  and  isotropic. 

The  following  assumptions  are  also  made  on  Equation  22. 

1.  steady-state, 

2.  no  internal  heat  generator,  and 

3.  the  temperature  variations  in  the  r 
and  0  directions  are  negligible. 

A  steady-state  assumption  is  made  since  the  material 
under  investigation  is  intended  for  use  as  an  insulator 
operating  under  steady-state  conditions. 

The  temperature  variations  in  the  r  and  9  directions 
are  considered  negligible  since  the  contact  region  is  small 
relative  to  the  particle  body. 

With  the  above  assumptions,  Equation  22  reduces  to  a 
one-dimensional  equation  and  heat  flow  through  a  contact 
region  is  viewed  as  quasi-one-dimensional.   Thus, 
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&  -  o  w 

A  method  of  solution  for  one -dimensional  problems  in 
which  q  =  constant  at  every  cross-section  is  [27] 

d[qxA(x)]  =  0  (24) 

which  gives 

QT  =  q  A(x)  =  constant  (25) 

The  x  denotes  heat  transfer  in  the  x-direction  and 

Introducting  Equation  26  into  Equation  25,  rearranging 
and  integrating  between  limits  gives 

TH-TT  T  -T 

Q  =  2 — k =  _g — b  (27) 

(l/ks)  /**  dx/A(x)      Ke 

where 

Re  =  (1/ks}  JxJ  dx/A^)  C28) 

and  represents  the  conductive  resistance  through  a  contact 
region.   Equation  28  shows  that  conductive  resistance  to 
heat  flow  is  directly  proportional  to  the  path  length,  dx, 
and  inversely  proportional  to  the  area,  A(x) . 

An  electron  microscope  photograph  of  a  diatomaceous 
earth  particle,  Figures  14  and  15,  page  92,  shows  that 
the  surface  of  these  particles,  where  contact  occurs,  has 
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rounded  or  smooth  geometric  features  as  shown  in  Figure 
4(a).   It  would  appear  that  the  regions  of  contact  could 
be  viewed  as  hemispherical,  cylindrical,  conical  and  so 
forth.   For  the  purpose  of  illustration,  a  hemispherical 
contact  region  will  be  assumed. 

A  contact  region  formed  by  two  hemispherical  pro- 
trusions is  shown  in  Figure  4(a).   Note  that  the  area  of 
contact  of  the  hemispherical  protrusions  has  been  flattened, 
as  illustrated  by  the  solid  line,  due  to  forces  imposed 
on  the  particles.   An  equivalent  cylindrical  heat  flow 
path  has  been  superimposed  on  the  hemispherical  contacts. 
This  is  to  illustrate  that  the  hemispherical  contact  region 
is  still  considered  quasi-one-dimensional.   The  use  of  a 
hemispherical  contact  region  serves  these  two  purposes 

1.  it  closely  resembles  the  actual  contact  region 
of  a  particle,  and 

2.  it  facilitates  the  calculations  of  the  contact 
area  with  the  use  of  the  Hertz  formula  [28]. 

The  symbols  used  in  analyzing  the  resistance  to  heat 
flow  through  a  contact  region  are  shown  in  Figure  4(b). 
Equation  28,  written  in  terms  shown  in  Figure  4(b)  is, 

Re  -  EJ  'x?  iyW  dx  C28) 

The  following  variable  transformations  are  made, 
y  =  R  sin  9 
x  =  R  cos  8 
dx  =  R  sin  9  d9 . 


Equivalent  cylindrical 
heat  flow  path 


particle  body 


(a) 

Parameters  Used  in  Analyzing  Heat  Flow  Through  A 
Hemispheric  Contact  Region 


Cb) 
Two  Hemispherical  Projections  Forming  A  Contact  Region 
Figure  4 
Two  Particles  Forming  A  Contact  Region 
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Substituting  into  Equation  28  and  simplifying, 


Re  ■  k^R  !l\  ike  «  C29) 

Evaluating, 

Re  =  FTR  ln  tan  2    I  A  (30) 

s  8 

and 

1  a  * 

Re  =  -  i n   In  tan  ^-  (31) 

k  ttR  2  ^  J 

s 

The  term  6   represents  the  angle  at  which  a  mimimum  contact 

radius,  y*,  occurs  for  a  given  particle.   The  contact  radius 

y*  is  defined  below.   Since  the  contact  angle  9*  is  less 

than  one  degree  as  calculated  using  the  Hertz  formula  (see 

Equation  35) ,  then  tan  8*/2  can  be  approximated  as 


tan  y  -  t  ~  y*/2(R-h)  (32) 
Now,  Equation  31  is 

Re  =  "  iAr  ln  TO^nT  (33) 
Since  R  >>  h,  the  conductive  resistance  can  be  written  as 

Re  =  i-^-n  In  ^  (34) 

k  ttR     y*  v  J 

s       ' 

The  value  of  y*  represents  the  radius  of  the  contact 

area  between  two  hemispheres  in  contact  and  it  is  calculated 
from  the  Hertz  formula  [28]  by 

y*  =  if    [Cl-H?)^  +  U  -  U?)/E2  ]}1/5               (35) 

"BT+  57 
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where 

F  =  load  per  contact 

Di  =  diameter  of  particle  one 

D2  =  diameter  of  particle  two 

y  =  Poisson's  ratio 

E  =  Young's  modulus 

In  Equation  35,  T>i    and  D2  are  diameters  of  spheres. 
In  this  text,  the  values  of  Di    and  D2  represent  the  dia- 
meters of  two  hemispherical  protrusions  that  form  a  contact 
point.   Hence,  the  parameters  T>i    and  D2  will  be  redefined 
as  D   and  D   ,  respectively.   Generally,  these  two 

cl        c2 

diameters  are  not  equal  but  throughout  this  text  the  values 

of  D   and  D   are  assumed  to  be  equal  and  will  therefore 

Ci       c2 

be  written  as  D  . 
c 

Equation  34  represents  the  conductive  resistance  to 
heat  flow  through  one  contact  region.   Combining  Equation 
34  with  Equation  27  and  rearranging  results  in  an  equation 
for  heat  flow  through  one  contact  region.   If  there  are 
Nt  particles  in  series  and  N.  particles  in  parallel,  then 
the  rate  of  heat  flow  through  one  foot  cube  of  material 
can  be  written 

n*  kg  Dc  NA(TH  -  Tc) 
QT  =  2Nh  ln(Dc/y*)  (36} 

where 
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n   =  number  of  contact  points  per  particle 
(heat  flow  in)  for  particle  in  series 

D   =  estimated  contact  region  diameter 

N.   =  number  of  particles  per  unit  of  area 
(parallel  heat  flow) 

N,   =  number  of  particles  in  a  unit  of 


'h 


height  (series  heat  flow) 


The  equivalent  contact  region  diameter,  Dc,  is  the 
estimated  diameter  of  a  hemispherical  protrusion  on  a  parti- 
cle.  This  was  shown  earlier  in  Figure  4(a).   Also  note  in 
Equation  36  that  radius,  R,  has  been  converted  into  terms 
of  diameter,  D  . 

For  a  unit  temperature  difference  across  a  given  depth 
of  insulating  material,  Equation  36  allows  one  to  define 
an  effective  coefficient  of  particle-to-particle  thermal 
conductivity  as 

im  k  D  N. 
v      -  s   c  A  r  37") 

Ks  "  2Nh  ln(Dc/y*J  ^/J 

The  above  equations  were  devised  on  the  basis  of  an 
array  of  uniform  particles.   However,  a  material  such  as 
diatomaceous  earth  consists  of  particles  that  are  non-uni- 
form in  size  and  shape  and  are  oriented  at  random.   There- 
fore, an  empirically  determined  correction  factor,  C,  must 
be  incorporated  into  Equation  37.   The  correction  factor 
will  take  into  account 

1.  bridging  of  particle  columns, 

2.  discontinuities  in  the  heat  flow  path,  such 
as  terminated  columns, 
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3.  non- uniformity  of  particle  size  and  shape,  and 

4.  changes  in  material  properties,  such  as  the 
addition  of  a  foreign  substance  in  a  pure 
base  material. 


In  view  of  the  above,  Equation  37  is  now 

C  7T  ks  n  Dc  NA 
Ks  =   2Nhln(Dc/y*)  (37) 

A  sample  evaluation  of  the  above  equation  will  be  discussed 
in  Chapter  V. 

Equation  37  was  derived  for  application  in  the  contact 
region  connecting  two  irregularly  shaped  particles.   As  a 
material  more  closely  approximates  a  sphere,  D  approaches 
the  diameter  of  the  sphere.   In  this  case,  N.  -  1/D  2  and 
N,  =  1/D  .   Equation  37  can  then  be  written  as 

C'k 
Ks  =  ln(Dc/y*)  (38) 

where 

C  =  Ctui/2 

The  particle-to-particle  conductivity,  K  ,  is  seen  to 
depend  primarily  on  the  particle  conductivity,  k  .   The 
natural  log  term  is  essentially  a  constant  for  typical 
variations  of  the  material  parameters  encountered  in  this 
study:  that  is,  for  a  given  powder  sample. 

The  value  of  C  is  a  constant  for  a  given  material. 
Incorporated  into  C  is  the  value  of  n,  the  number  of 
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contact  regions  leading  into  a  particle.    The  value 

of  K   is  seen  to  vary  directly  with  n.   Thus,  one  should 

s 

choose  an  insulating  powder  so  that  n  is  minimized. 


Radiant  Transfer  in  a  Porous  Medium 

Introduction.   The  ratio  of  radiant-to-total  heat 
transfer  in  a  porous  medium  has  been  modeled  extensively. 
However,  there  is  little  explicit  experimental  data, 
identifying  radiation  separately,  with  which  to  judge 
the  effectiveness  of  these  models.   In  a  porous  medium 
in  which  the  particles  are  opaque  and  small  (<10~2  ft) 
and  in  which  the  temperatures  are  less  than  2700°F,  there 
is  agreement  in  the  literature  that  radiant  compared  to 
total  heat  flow  is  negligible  [29,30]. 

Radiant  heat  transfer  in  a  porous  powder.   In  view  of 
the  above  paragraph,  net  radiant  transfer  through  powders 
for  conditions  encountered  in  this  study  was  considered 
to  be  negligible.   The  powders  used  were  on  the  order  of 
lO'4  ft  in  diameter  and  upper  temperature  limit  was  1000°F. 
As  will  be  demonstrated  in  Chapter  V,  net  radiant  transfer 
through  diatomaceous  earth  powder  under  typical  operating 
conditions  accounted  for  approximately  1% ,  on  a  theoretical 
basis,  of  total  energy  transfer.   Even  though  this  is 


Contact  regions  for  energy  transport  out  of  a  particle  are 
considered  as  a  contact  region  leading  into  an  adjacent 
particle. 
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considered  to  be  negligible,  a  theoretical  analysis  of 
radiant  heat  transfer  through  powders  will  be  presented 
for  completeness. 

The  theoretical  model  presented  is  similar  to  those 
of  Argo  and  Smith  [29]  and  Schotte  [30].   To  model  radiant 
transport  the  above  authors  assumed  spherical  particles 
and  considered  radiant  transfer  through  parallel  planes 
located  on  each  side  of  a  particle.   See  Figure  5. 


6    Ti  >  T- 


Figure  5 

Model  for  Derivation  of  Radiant  Contribution  to 
Thermal  Conductivity  through  a  Porous  Medium 


The  model  considers  radiant  transfer  through  the 
voids  between  particles  and  radiant  transfer  from  surface 
to  surface  in  series  with  the  solid  particle  conduction. 


46 


The  model  also  assumes  that  the  particles  were  opaque 
and  that  the  particle  size  was  large  compared  to  the 
wavelength  of  radiation.   Considering  the  last  restriction, 
the  particle  size  used  in  this  study  was  15  -  125  ym. 
For  a  temperature  range  of  250  -  1000°F,  Wien's  displace- 
ment law  gives  an  average  wavelength  of  7.35  -  3.57  urn, 
respectively.   Thus,  the  condition  is  well  satisfied 
with  large  particles  at  all  conditions,  reasonably 
well  satisfied  with  medium  size  particles  at  higher 
temperatures  or  small  particles  at  lower  temperatures, 
but  is  doubtful  for  the  smallest  particles  at  the  lowest 
temperatures . 

One  may  view  radiant  energy  as  a  packet  of  energy, 
or  photon,  instead  of  a  wave.   Of  the  many  photons  that 
participate  in  radiation  passing  through  a  porous  powder 
most  will  strike  the  surfaces  that  form  a  void  many  times 
before  they  could  escape.   As  a  result  the  probability  of 
absorption  is  very  high  compared  to  the  probability  of 
escape.   Hence,  the  effective  absorptivity  and  emmissivity 
may  be  approximated  as  unity.   By  viewing  radiant  transfer 
in  terms  of  photons,  it  is  easy  to  visualize  radiant  energy 
transfer  as  meeting  great  resistance  in  flowing  through  a 
porous  powder,  especially  for  conditions  previously 
described. 

The  results  of  the  analysis  by  Schotte  [30]  give  an 
overall  radiative  conductivity  coefficient  of 
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Kr  -  T-i-^i   ♦  akrv  (29) 

*7    ^ 

where  6  is  the  void  fraction,  kg  is  the  solid  particle 
conductivity,  and  k   is  a  particle  radiation  conductivity 
coefficient.   The  latter  variable  represents  radiant  heat 
transfer  between  a  particle  and  its  neighbor.   This 
value  will  be  defined  shortly. 

The  right-hand  side  of  Equation  39  consists  of  two 
terms.   The  first  term  represents  radiant  heat  transfer  to 
a  particle  in  series  with  conduction  through  the  particle. 
The  second  term  accounts  for  radiation  through  the  void 
space  adjacent  to  the  particle. 

The  radiative   coefficient,  kry,  which  represents 
radiative  transfer  between  a  particle  and  its  neighbor, 
is  given  by  Schotte  [30]  as , + 

krv  =  4ed?T3  (40) 

where  e  is  the  emissivity  of  the  particles,  3"  is  the 
average  void  dimension,  £  is  the  Stef an-Boltzmann  con- 
stant, and  T  is  absolute  temperature.   The  temperature 
used  in  this  equation,  T,  is  the  average  temperature  of 
the  particles.   Argo  and  Smith  [29]  have  pointed  out  that 


+  This  equation  is  presented  by  Schotte  [30]  as  kry  ■  0.692d 
T3/108.   This  equation  and  the  one  used  above  are  the  same. 
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one  investigator  [31]  used  a  bulk  mean  temperature  in 
evaluating  Equation  36  and  found  it  in  agreement  with  more 
elaborate  methods  for  determining  radiative  transfer 
through  powders. 

In  summary,  the  above  paragraphs  illustrate  the 
role  of  radiant  to  total  energy  transport  through  a  porous 
powder.   The  parameters  that  apparently  control  radiant 
heat  flow  through  powders  are  the  particle  size,  opacity 
of  the  material,  particle  conductivity,  void  dimension, 
and  boundary  temperatures. 


CHAPTER  IV 

EQUIPMENT  DESIGN 
AND  MATERIAL  PREPARATION 


Introduction 

The  intent  of  the  experimental  procedure  was  two- 
fold.  First,  apparent  thermal  conductivities  for  various 
particle  sizes  were  needed  to  establish  the  required 
insulation  thicknesses  for  a  high  temperature  storage 
system.   This  required  the  design  and  construction  of 
a  conductivity  test  apparatus  capable  of  measuring  con- 
ductivity of  a  powder  material  with  a  hot  face  temperature 
up  to  1000°F  and  for  vacuum  levels  down  to  100  ym  Hg. 

The  second  objective  of  the  experimental  procedure 
was  to  design  a  prototype  high  temperature  thermal  energy 
storage  system  using  a  porous  powder  under  a  moderate 
vacuum  (P  ~  1  mm  Hg)  as  the  insulating  material.   The 
thermal  energy  storage  vessel  was  designed  to  use  an 
organic  oil  (corn,  cottonseed,  peanut,  or  soybean)  as  the 
storage  medium. 

Conductivity  Test  Apparatus 

The  conductivity  test  apparatus  is  shown  in  Figures 
6(a)  and  6(b).   The  central  core  of  the  test  cell  consists 
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electrical  wiring  port 

wax  cup  thermocouple  port 
V  OD  pipe 


acuum  port 


centering  plate 


2-5-  OD  copper  sleeve 


vacuum  gauge 
port 


1~-  aluminum  plug, 
four  threads 
per  inch 


1" 
4-q   OD  copper 

cylinder 


centering  plate 
vacuum  port 


steel  wool  (grade  000) 
(a) 
Sectional  View  Of  Conductivity  Test  Apparatus 

Figure  6 
Conductivity  Test  Apparatus 
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Top  View 


Centering  Plate 
Figure  6(a)  (Cont'd.) 
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(b) 
Figure  6 


Photograph  of  Conductivity  Test  Apparatus 
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of  two  guard  heaters  and  a  center  test  section  heater. 
The  nichrome  heater  wires  were  wound  on  a  threaded 
aluminum  core  which  served  to  thermally  dampen  small  devia- 
tions in  the  power  supply.   The  aluminum  core  was  wrapped 
with  a  ceramic  paper  and  inserted  into  a  copper  sleeve. 
The  copper  sleeve  helped  to  insure  a  uniformly  heated 
surface.   Chromel-alumel  thermocouples  were  welded  to  the 
copper  sleeve  so  that  a  temperature  difference  between 
the  guard  heaters  and  the  center  test  section  heater 
could  be  read.   In  this  manner,  the  guard  heaters  could 
be  adjusted  to  the  same  temperature  (+  1/2°F)  as  the 
center  test  section  heater,  insuring  radial  heat  flow  in 
that  section.   By  measuring  the  electrical  energy  dissi- 
pated in  the  central  test  section,  and  by  measuring  the 
surface  temperatures  at  the  boundaries  of  the  powder 
insulation,  the  apparent  thermal  conductivity  was  calcula- 
ted.  For  a  cylindrical  geometry, 


2irK    (TH-Tf)L 

Q  =     _a H L (4-g 

xelec  lnfr   /r-J 

Solving  for  conductivity, 

Q  ,    ln(r  /r.) 
k     -   xelec    ^  o'  3/  r/m 

Ka  "   2tt  (Th-Tc)L  t4ZJ 


For  this  apparatus, 
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Qelec   =   rate  of  electrical  energy  dissipation 

L         1  ft  test  section  length 

rQ     =   0.163  ft  outside  radius  of  insulating 
material 

ri     =   0.089  ft  inside  radius  of  insulating 
material 

Tc     =   average  cold  surface  temperature  of  insu- 
lation at  boundary 

TH     =   average  hot  surface  temperature  of  insu- 
lation at  boundary 

The  electrical  energy  dissipated  was  calculated  by  measuring 
the  voltage  and  current  of  the  center  test  section  heater 
wires.   Thus, 

Qelec  =  EI  cos  *  (43) 

where, 

cos  0=1  (assuming  the  power  factor  to  be  unity 
for  resistance  heating) 

E     =  potential 

I     =  current 

The  value  of  Tc,  the  cold  wall  temperature,  was  measured 
on  the  outside  surface  of  the  conductivity  test  cell. 
The  calculated  estimate  of  the  temperature  drop  across 
the  copper  wall  was  negligible  (<  0.002°F).   The  heat 
losses  through  the  wires,  thermocouples, and  electrical 
conductors  had  a  negligible  effect  on  the  test  section 
energy  balance,  since  the  temperature  difference  between 
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that  section  and  the  guard  heaters  was  held  to  +  1/2°F 
of  each  other. 

Pads  of  fine  steel  wool  were  placed  at  the  top  and 
bottom  of  the  conductivity  test  cell  to  facilitate  evacua- 
ting the  fine  powder  without  entrainment .   The  vacuum 
ports  at  the  top  and  bottom  of  the  test  cell  permitted 
a  faster  evacuation  time  than  a  single  port. 

A  vacuum  port  was  placed  above  the  upper  guard 
heater  and  below  the  lower  guard  heater  so  that  its 
presence  would  not  interfere  with  the  center  test  section. 
Several  different  types  of  vacuum  gauges  were  used  during 
the  course  of  testing. 

The  first  was  a  thermocouple  gauge  with  a  range  of  one 
atmosphere  to  1  ym  Hg.   However,  it  was  difficult  to 
maintain  calibration  due  to  powders  infiltrating  its 
sensing  element  and,  consequently,  this  gauge  was  not  used. 
A  swivel  type  McLeod  gauge  with  a  range  of  1  ym  -  5  mm  Hg 
was  the  second  type  used.   The  gauge  worked  well  after  a 
small  piece  of  fine  steel  wool  (0000)  was  put  in  line  with 
it  to  prevent  powders  from  reaching  the  mercury. 

Alignment  of  the  central  core  of  the  test  apparatus 
was  maintained  by  a  centering  plate  as  shown  in  Figure  6(a) 

Wire  feeds  were  created  by  packing  asbestos  string 
around  the  wires.  The  wires  were  sealed  as  they  passed 
into  the  test  apparatus  by  melting  high  temperature  wax 
around  them.   The  wax  had  to  be  re-heated  periodically 
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due  to  differential  expansion  between  the  wire  and 


The  Design  of  a  High  Temperature  Thermal  Energy  Storage 
System 

A  high  temperature  thermal  energy  storage  system  was 

designed  according  to  the  following  criteria. 

1.  Heat  loss  per  day  not  to  exceel  10  Btu/hrft2. 

2.  The  inner  tank  has  a  minimum  high  pressure 
rating  of  150  psi  at  700°F. 

3.  Both  the  inner  and  outer  tank  have  low  pressure 
rating  at  1.32  x  10-5  psi  (100  ym  Hg) . 

4.  The  inner  tank  is  capable  of  holding  organic 
oils . 

5.  The  vessel  construction  is  simple  and  durable. 

6.  The  vessel  requires  little  maintenance. 

7.  The  inner  storage  vessel  is  sized  to  be  compati- 
ble with  a  6'  x  8'  concentrating  collector. 

8.  The  design  potential  is  intended  for  mass  pro- 
duction at  an  acceptable  cost. 

The  tank  design  is  shown  in  Figure  7. 

The  inner  tank  was  made  from  a  14  gauge  galvanized 
steel,  42  gallon  water  tank.   Note  that  one  end  of  the 
water  tank  was  changed  from  concave  to  convex  to  insure 
stability  at  elevated  pressures. 

The  outer  tank  was  constructed  from  a  14  gauge 
galvanized  steel  200  gallon  water  tank.   The  cylindrical 
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5/16"  bolts,  typical 


Figure  7(b)  (Cont'd) 

End  View  of  Tank  Showing  Internal 
Cable  Support  and  Flange  Bolt  Pattern 
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1" 
—  x  2' 

16  x  ; 

angle  iron 


x  2' 


J  x  1"  x  1" 
angle  iron 


Figure  7(c)  (Cont'd) 
External  Support  for  Inner  Tank 
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(d) 

Figure  7 

Photograph  of  Prototype  High  Temperature 
Thermal  Energy  Storage  System 


61 

portion  of  the  tank  was  shortened  and  the  concave  end  was 
changed  to  convex  so  that  the  insulation  thickness  at 
the  ends  of  the  tanks  would  be  equal  to  the  thickness 
between  the  cylindrical  portions  of  the  tank. 

The  inner  tank  was  suspended  by  two  1/16"  diameter 
stainless  steel  cables.   These  cables  were  selected  for 
their  high  breaking  strength  (500  lb.  at  room  temperature) 
and  low  thermal  conductivity  (10  Btu/hrft°F  at  room  tempera- 
ture) .   The  two  cables  had  sufficient  strength  at  elevated 
temperatures  (1000°F)  to  afford  a  safety  factor  of  two  (2). 

The  inlet  and  outlet  pipes  were  designed  with  a 
U-shaped  trap  which  served  as  an  expansion  loop  and  a 
heat  trap.   The  inlet  pipe  at  the  top  had  the  loop  placed 
close  to  the  inner  tank  so  that  it  would  prevent  volatile 
vapors  from  the  organic  oil  from  entering  the  outlet  pipe 
when  the  pipe  was  not  in  use.   The  loop  on  the  botton 
outlet  pipe  was  placed  at  the  opposite  end  of  its  point 
of  attachment  and  the  outlet  pipe  was  sloped  downward  so 
that  heat  conducted  through  the  metal  would  induce  convective 
flow  in  the  opposite  direction,  thereby  helping  to  minimize 
heat  flow  out  of  the  inner  tank. 

The  rings  on  the  outer  tank  were  needed  to  meet 
safety  recommendations  on  externally  pressurized  tank 
design  [32]. 

The  inner  tank  was  designed  to  contain  organic  oils. 
Copper  metal  cannot  be  used  with  organic  oils  since  these 
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oils  are  acidic  in  nature  and  react  chemically  with  copper. 
The  long  term  effect  of  high  temperature  organic  oils  on 
steel  is  presently  unknown.   At  low  temperatures,  steel 
containers  have  been  used  successfully  to  contain  organic 
oils . 

The  minimum  storage  vessel  size  was  estimated  for 
a  6  ft  x  8  ft  (aperture  area)  concentrating  collector 
as  follows. 

Heat  delivered  to  _  Heat  stored  in  oil  +  Heat  losses     ,..>. 
storage  system      and  its  container    from  piping    ^  * 

Assuming  the  pipe  losses  equal  10%  of  the  energy  stored 
in  the  oil  and  its  container,  then 


H.   =  1.10  H  ..  (45) 

in         sto  v  J 


and 


IsAantx  =  l.lpVCpAT  (46) 


where  I   is  solar  insolation,  A  is  the  collector  aperture 
area,  n   is  the  collector  average  thermal  efficiency,  p  is 
density  of  organic  oil  (cottonseed) ,  V  is  the  oil  tank 
volume,  C   is  the  average  specific  heat  of  the  oil,  and 
AT  is  the  temperature  rise  of  the  oil. 
Solving  for  volume 


I  A  n  x 

V  =   s  a  z T471 

1.1  C  AT  ^/} 
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Assuming  the  following  values, 
n   =   50% 

T 

Is  =   300  Btu/hrft2 

x  =8  hr/day 

P  =   46.5  lbm/ft3  @  600°F  (estimated)* 

C  =   0.66  Btu/lbm°F  @  600°F  (estimated)* 

AT  =   400°F  (200  -  600°F) 

A   =   48  ft2 
a 

and  solving  for  V  in  Equation  47  yields 

V   =   4.27  ft3  (maximum) 

A  42  gallon  tank  holds  5.62  ft3  and  afforded  an  excess 
capacity  to  serve  as  an  expansion  chamber. 

The  outer  tank  size  was  chosen  from  commercially 
available  tanks.   The  most  convenient  size  was  a  220 
gallon  galvanized  water  tank.   The  outer  tank  diameter 
was  2.5  ft  and  the  inner  tank  outside  diameter  was  1.33 
ft.   Thus,  there  was  a  0.585  ft  space  for  insulation.   The 
ends  of  the  tank  were  designed  to  leave  an  equivalent  0.585 
ft  space. 

A  preliminary  test  on  diatomaceous  earth  powder  indi- 
cated that  an  apparent  conductivity  value  of  K   =0.01 

a 


Cottonseed  oil  data  [3] 
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Btu/hrft°F  could  be  obtained  at  a  hot  surface  temperature 
of  600°F,  a  cold  surface  temperature  of  70°F,  and  at  1  mm 
Hg  pressure.   A  0.585  ft  thick  layer  of  insulating  powder 
for  these  conditions  would  hold  the  daily  energy  loss 
of  the  storage  tank  to  approximately  9%  per  day  of  the 
initial  energy  stored.   This  approximation  is  based  on 
changes  in  temperature  and  average  values  of  C  over  the 
interval  of  temperature  change. 

Diatomaceous  Earth  Properties 

Diatomaceous  earth  is  the  fossilized  remains  of  micro- 
scopic marine  algae.   It  consists  of  approximately  98% 
Si02  (silica)  .   The  diatomaceous  earth  used  in  this  study 
was  purchased  commercially.   Its  intended  purpose  was  for 
use  in  swimming  pool  filters.   This  material  has  been 
used  industrially  as  high  temperature,  loose  fill,  insula- 
tion and  in  compressed  brick  form  as  a  furnace  lining. 
The  loose  fill  material  has  an  upper  temperature  limit  of 
1600°F  and  the  brick  has  an  upper  limit  of  2500°F  [19]. 
The  apparent  reason  for  the  limitation  of  the  loose  fill 
material  is  that  the  powder  particles  begin  to  soften 
at  that  temperature.   This  would  increase  the  point-to- 
point  solid  conductivity  and  also  cause  loss  of  volume  of 
the  material  due  to  settling. 

In  fine  powder  insulation  it  is  useful  to  determine 
the  effect  of  particle  size  on  apparent  thermal  conductivity, 
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One  difficulty  of  determining  particle  size  is  that  fine 
powders  tend  to  clump.   An  explanation  for  this  phenomenon 
will  be  discussed  later.   However,  this  made  it  difficult 
to  separate  these  particles  in  a  sieve  and  shaker.   Grind- 
ing these  particles  in  a  ball  and  mill  reduced  their  size 
but  did  not  make  it  easier  to  separate  them  into  various 
particle  sizes. 

It  was  found  that  the  addition  of  small  amount  of 
fine  carbon  black  to  the  diatomaceous  earth,  either  during 
the  separation  process  or  the  grinding  process,  permitted 
the  powder  to  separate  quickly  in  a  sieve  and  shaker.   The 
reason  for  this  is  that  the  fine  carbon  particles  act  like 
a  lubricant.   Carbon  particles  are  shown  in  Figure  11. 
Note  the  puffy  appearance  of  the  agglomerated  carbon 
particles.   This  suggests  that  the  large  particles  may  be 
comprised  of  many  smaller  carbon  particles. 

A  possible  explanation  of  how  carbon  acts  as  a  lubri- 
cant may  derive  from  the  fact  that  carbon  is  electrically 
conductive.   When  the  diatomaceous  earth  particles  are 
vibrated,  the  particles  become  charged  with  static  electri- 
city.  Hence,  if  particles  have  opposing  charges,  they 
attract  each  other.   This  prevents  the  particles  from  pass- 
ing through  the  sieve  and,  consequently,  no  separation  of 
particles  occur.   When  carbon  powder  is  added  to  diato- 
maceous earth,  it  short  circuits  the  static  charge  between 
the  particles  and  this  permits  the  particles  to  pass  through 
the  sieve. 
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After  the  diatomaceous  earth  was  separated  into  its 
various  particle  sizes,  its  geometry  remained  essentially 
the  same.   Figure  8  shows  particles  trapped  in  Tyler 
Sieve  #500.   Figure  9  shows  particles  trapped  by  Tyler 
Sieve  #170.   Both  of  these  photographs  vividly  display 
the  irregular  characteristics  of  these  powders. 

The  various  particle  size  distributions  obtained  in 
a  sieve  and  shaker  are  given  in  Table  3.   The  Tyler  sieve 
number  corresponds  to  screen  weaves  of  the  following  opening 
size 

Tyler  Sieve  No.     120   170   200   270   325   500   <500 
Opening  (microns)   125   90    75    53   45    25   <  25 

The  last  category,  <500,  was  simply  a  catch  pan  for  particles 
falling  through  the  #500  screen. 

The  values  shown  in  Table  3  are  representative  only 
of  the  particle  sizes  in  the  original  diatomaceous  earth. 
These  values  changed  somewhat  with  each  test. 

The  purpose  of  generating  data  for  Table  3  was  to 
find  the  various  nominal  particle  sizes  with  their  corres- 
ponding quantities.   The  variables  influencing  the  separa- 
tion process  were  the  length  of  time  the  shaker  was  operated 
and  the  amount  of  carbon  present.   The  objective  was  to 
use  as  little  carbon  as  possible.   The  general  results 
of  Table  3  follow. 
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Electron  Microscope  Photograph  of  Diatomaceous 

Earth-Carbon  (II  wt)  Mixture,  Particles  from 

Tyler  Sieve  #5  00  (x  1,000 


Figure  9 


Electron  Microscope  Photograph  of  Diatomaceous 
Earth-Carbon  (II  wt)  Mixture,  Particles  from 
Tyler  Sieve  #170  (x  1,000) 
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Figure  10 

Electron  Microscope  Photograph  of  Diatomaceous  Earth 
(as  purchased),  x  2,000 


Figure  11 

Electron  Microscope  Photograph  of  Carbon  Black, 
Monarch  500  (x  10,000) 
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1.  Run  #1  illustrates  the  difficulty  of 
separating  raw*  diatomaceous  earth. 

Even  after  30  minutes,  very  few  particles 
were  separated. 

2.  Run  #2  illustrates  that  a  large  percentage 
of  the  particles  in  the  powder  are  indeed 
small . 

3.  The  general  trend  is  that  as  the  amount  of 
carbon  increased,  the  greater  the  percentage 
of  small  particles  and  the  faster  they 
separated. 

4.  A  given  type  of  powder  with  the  same  amount 
of  carbon  and  the  same  shake  time  shows  a 
general  consistency  in  particle  distribution, 
but  some  particles  of  a  given  size  may  vary 
appreciably.   Compare  Run  #4  and  Run  #5. 

5.  Carbon  added  prior  to  pulverzing  diatoma- 
ceous earth  permits  the  separation  of 
fine  particles  in  less  time  than  when  it 
was  added  after  pulverizing. 


Representative  data  on  density  of  the  various  particle 

sizes  are  given  in  Table  4.   Note  the  discrepancy  between 

f 
the  measured  void  fractions   and  the  calculated  void 

fractions.   This  could  be  due  to  air  being  trapped  in  and 

among  the  particles  and  not  being  displaced  during  the 

void  measurement  process.   Another  explanation,  in  view  of 

Figures  8  and  9,  is  that  the  particles  do  not  separate  well 

and  that  the  mixture  of  large  and  small  particles  causes  the 

void  fraction  to  remain  relatively  constant. 


A 

Raw  diatomaceous  earth  refers  to  the  material  as  purchased, 

Void  fraction  is  the  ratio  of  the  void  space  in  a  con- 
tainer to  the  total  volume  of  the  container. 
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The  method  used  to  measure  the  void  fraction  and 
density  was  to  fill  carefully  a  graduated  cylinder  of 
known  weight  with  the  powder  and  gently  tap  the  cylinder 
to  settle  the  particles.   The  cylinder  with  the  particles 
was  then  weighted.   The  cylinder  with  the  powder  was  then 
filled  with  water  to  the  original  level  of  the  powders, 
and  it  was  then  weighed  again.   With  a  knowledge  of  the 
volume  and  weight  of  the  powders,  a  density  was  then  cal- 
culated.  The  void  fraction  was  found  by  deducing  the  volume 
of  water  required  to  fill  the  graduated  cylinder  to  the 
original  volume  of  powder.   The  ratio  of  the  latter  two 
quantities  is  the  void  fraction. 


CHAPTER  V 
EXPERIMENTAL  PROCEDURE  AND  RESULTS 

Conductivity  Test  Equipment 

The  conductivity  test  equipment  arrangement  is  shown 
in  Figures  12  and  13.   Power  was  supplied  through  rheostats 
connected  to  a  standard  120  volt  AC  outlet.   Variations  in 
the  AC  outlet  due  to  daytime  peaks  were  overcome  by  obtain- 
ing a  long-term  average  for  the  data  over  a  period  of 
several  days  to  more  than  a  week  for  each  datum  point. 
These  data  were  recorded  manually. 

The  test  section  heaters  were  monitored  for  current 
and  voltage.   This  was  accomplished  by  placing  an  ammeter  on 
one  lead  of  the  heater  wire.   Voltage  was  measured  on  the 
heater  side  of  the  power  supply.   The  guard  heater  power 
was  monitored  only  for  interest  in  their  power  levels 
relative  to  the  central  heater.   The  guard  heater  power 
levels  were  adjusted  to  give  a  temperature  within  +0.5°F 
of  the  central  heater.   This  power  level  was  usually  slightly 
more  than  the  power  level  of  the  central  heater. 

The  temperature  of  the  four  thermocouples  inside  the 
test  cell  were  measured  with  a  digital  volt-ohm  meter. 
An  ice  bath  was  used  as  a  temperature  reference  for  the 
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Figure  12 
Conductivity  Test  Equipment  Arrangement 
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Figure  13 
Conductivity  Test  Equipment 
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chromel-alumel  (Type  K)  thermocouples.   The  temperature 
levels  were  recorded  manually. 

A  mercury-filled  McLeod  gauge  with  a  range  of  1  urn  -  5 
mm  Hg  was  used  to  measure  the  vacuum  level  in  the  conducti- 
vity test  cell  for  data  presented  in  this  study.   For  the 
preliminary  testing  period  a  thermocouple  vacuum  gauge 
was  used.   It  was  found  to  lose  calibration  often  and  a 
standard  gauge  was  needed  for  recalibration.   It  was  found 
easier  to  use  the  standard  gauge  (McLeod)  directly  in  the 
experiment . 

A  common  vacuum  pump  with  a  blank-off  pressure  of  15 
ym  Hg  was  used  for  evacuating  the  system.   A  variable 
leak  valve  was  placed  at  the  bottom  of  the  test  cell  to 
regulate  pressure. 

The  objective  of  this  experiment  was  to  find  the 
change  in  thermal  conductivity  due  to  changes  in 

1.  particle  size, 

2.  vacuum  level,  and 

3.  temperature  levels. 

Thermal  Conductivity  of  Powder  Insulation  Under  Vacuum 

For  the  determination  of  the  thermal  conductivity  of 
a  powder,  the  thermal  conductivity  apparatus,  as  shown  in 
Figure  6,  was  filled  with  the  powder  to  be  tested.   The 
powders  were  heated  prior  to  evacuation  so  that  the  elevated 
temperatures  would  aid  moisture  removal  from  the  powder 
particles . 
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Once  the  system  was  thoroughly  heated,  it  was  sealed 
and  evacuation  was  initiated.   The  heater  elements  were 
adjusted  so  that  the  main  heater  and  guard  heater  thermo- 
couples read  within  1/2°F  of  each  other.   This  process 
usually  required  several  days.   After  a  uniform  temperature 
was  reached,  the  system  was  generally  operated  under  this 
condition  for  at  least  one  day  to  insure  that  steady-state 
had  been  achieved.   This  additional  time  was  necessary 
to  minimize  the  effect  of  the  daily  changes  in  the  AC 
power  supply  and  the  daily  temperature  swing  of  the  test 
room. 

During  the  testing  procedure,  data  were  taken  manually, 
usually  on  hourly  intervals,  for  the  three  heater  tempera- 
tures (6  points),  outside  surface  of  test  apparatus,  test 
heater  potential  and  current,  and  room  temperature. 

Data  and  Results  from  Thermal  Conductivity  Test 

The  experimental  values  obtained  from  the  thermal 
conductivity  apparatus  are  shown  in  Table  5.  The  four 
measurements  were 

1.  pressure  level,  measured  in  ym  Hg, 

2.  the  hot  surface  temperature,  Th,  bounding  the 
powder  material  adjacent  to  the  main  heater, 

3.  the  cold  surface  temperature,  Trj,  taken  on 
the  outside  surface  of  the  test  apparatus, 

4.  power  dissipated  from  the  central  test  section 
heater  (potential  x  current). 
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From  the  measured  data  the  apparent  thermal  conductivity 
coefficient,  K  ,  was  calculated  with  the  use  of  Equation 

42. 


Discussion  of  Results  of  Apparent  Conductivity  Test 

Range  of  values.   The  range  of  apparent  conductivity 
values,  as  seen  in  Table  5,  vary  as  follows 


Highest  value:   K  =  0.0648  Btu/hrft°F 

3- 

@  TH  =  631°F,  P  =  760  mm  Hg 


Lowest  value:    K   =  0.0029  Btu/hrft°F 

a  ' 

@  TH  =  340°F,  P  =  0.6  mm  Hg 


To  put  these  values  in  perspective,  air  has  the  following 
experimental  conductivity  values  [33] 

Air  at  631°F:    K  =  0.0272  Btu/hrft°F 
@  P  =  760  mm  Hg 

Air  at  340°F:    K  =  0.0212  Btu/hr  ft°F 
@  P  =  760  mm  Hg 

From  this,  the  highest  test  data  value,  taken  at  760  mm  Hg 
pressure,  is  2.4  times  higher  than  the  conductivity  of  air, 
but  the  lowest  value,  with  a  pressure  of  0.6  mm  Hg,  is  7.3 
times  less  than  air  at  760  mm  Hg.   This  is  significant  since 
the  conductivity  of  air  is  essentially  independent  of 
pressure  except  under  the  low  pressure  conditions  (Kn  >  1) 
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TABLE  5 

Average  Thermal  Conductivity  Values 
For  Evacuated  Diatomaceous  Earth* 


Powder  in  Tyler  Sieve  #500  (loose  powder) 


Pressure 
(urn  Hg) 

TH 
(°F) 

TC 
(°F) 

Ka 
(Btu/hrft°F) 

1 

500 

405 

90 

.00364 

2 

600 

733 

101 

.00648 

3 

800 

700 

93 

.00750 

4 

850 

330 

89 

.00483 

5 

1000 

730 

104 

.00722 

6 

1000 

660 

85 

.00790 

7 

1000 

558 

97 

.00444 

8 

1000 

452 

87 

.00407 

9 

1000 

443 

84 

.00417 

10 

760,000 

631 

192 

.06480 

11 

760,000 

432 

154 

.04600 

Powder  in  Tyler  Sieve  #500  (powders  physically  compressed 
to  45  psi) 


12 

600 

13 

725 

14 

650 

15 

725 

16 

775 

17 

800 

18 

850 

19 

900 

20 

1000 

21 

1000 

340 

97 

.00290 

472 

85 

.00380 

498 

100 

.00540 

691 

100 

.00910 

295 

86 

.00338 

300 

84 

.00327 

685 

103 

.00910 

677 

102 

.00930 

694 

108 

.01230 

746 

111 

.01130 

"JT 

Powders  contain  approximately  1%    carbon  by  weight, 
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TABLE  5,  Continued 


Powder  in  #325  pan  (loose  powder) 

Pressure      Th  Tc  Ka 

(pm  Hg)       (°F)  (°F)      (Btu/hrft  F) 

22  1400  740  100        .00580 

23  1800  462  90         .00475 

Powder  in  #170  pan  (loose  powder) 


24 

175 

25 

200 

26 

300 

27 

800 

28 

1000 

29 

1000 

30 

1000 

31 

3800 

456 

81 

.00635 

699 

88 

.00656 

588 

108 

.00590 

647 

110 

.00723 

344 

80 

.00972 

531 

8  5 

.01198 

710 

98 

.01343 

557 

95 

.01714 

Raw  Diatomaceous  Earth  (No  carbon  added  to  loose  powder) 


32 

600 

33 

1000 

34 

2000 

35 

8000 

36 

8000 

441 

93 

.00920 

725 

110 

.01200 

545 

96 

.00980 

386 

110 

.02900 

583 

117 

.02740 

Pulverized  Diatomaceous  Earth  (loose  powder) 


37 

1000 

38 

8000 

39 

8000 

Air  Only 

40  760,000 
Fiberglass  Cloth 

41  400 


446 

99 

.00960 

511 

101 

.01560 

715 

116 

.01448 

603 

235 

.07800 

inum  Foil 

■  5  lay 

ers 

720 

110 

.01440 

discussed  earlier.   This  indicates  that  a  powder  under 
vacuum  is  effective  insulation  at  elevated  temperatures. 

The  particular  area  of  interest  in  this  study  was  the 
conductivity  values  of  the  powders  near  1  mm  Hg,  since 
this  vacuum  level  was  within  the  range  of  common  mechanical 
vacuum  pumps.   The  conductivities  of  several  of  the  more 
thoroughly  investigated  powders  are  shown  in  Figure  16 
(page  96).   This  graph  shows  the  approximate  limits  and 
the  general  trend  of  the  apparent  conductivities  associated 
with  the  large  and  small  grain  powders  tested. 

Radiative  transfer.    The  net  contribution  due  to 
radiation  through  a  powder  insulation,  in  the  temperature 
range  of  interest  in  this  study,  appears  to  be  negligible. 
Recall  Equations  39  and  40, 


Kr  -  l/ksVl/krv  +  *  *rv  ^ 


k    =  4e  d~  ?  T3  (40) 

rv 


Using  typical  values  in  these  equations,  the  resulting 
coefficient  of  radiative  conductivity  is  small  compared 
to  other  modes  of  energy  transfer.   For  example,  using  the 
following  values  for  particles  in  Tyler  sieve  #500 
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5  =0.76  (Table  4) 

e  =1.0  (assumed  value) 

d  -  0.5D  =  5.75  x  10~D  ft. 

k  =0.76  Btu/hrft°F  (pyrex  glass) 

and  calculating  for  radiant  transfer  for  powder  in  a 
plane  adjacent  to  a  hot  surface  boundary  of  1000 °F 
(1460°R)  and  for  a  plane  at  250°R  gives  a  range  of 


K   =  0.001  to  0.0001  3tu/hrft°F 
r 


Comparing  the  above  values  to  a  typical  apparent  conductivi- 
ties of  K   =  0.015  and  0.005  Btu/hrft°F,  shows  that  radiation 
a  ' 

accounts  for  7%  and  2%  of  the  totals,  respectively.   These 
values  represent  a  worst  case  for  both  temperatures  since 
they  were  calculated  for  radiant  transfer  near  the  hot 
surface  boundary.   For  powder  near  a  cold  face  boundary  at 
T   =  90°F  (550°R) ,  then  K^  represents  about   1%   compared 
to  an  average  value  of  K   =0.01  Btu/hrft°F.   Thus,  the 
net  contribution  of  radiant  transfer  through  a  powder  for 
the  temperature  limits  of  this  study  appears  to  be  negligible. 

Effective  gas  conductivity.   Note  that  Equation  20, 
when  rearranged  gives 
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(V5)   2Ru  h       \  /2  *    °2  AN  -1 

Kg  =  4X7^17  W  (^  +  iTT } 

Evaluating  with  the  following  parameters  for  air  alone  in 
a  container  of  one  foot  cubed, 


p 

= 

760,000  m  Hg 

T 

= 

492°R  (32°F) 

M 

= 

28.97  lb  /lb  -mole 
m   m 

a 

= 

1.1975  x  10"9  ft 

d 

= 

1.0  ft  (since  d"  >>  A) 

Y 

= 

1.39 

c„ 

= 

0.18  Btu/lb  °F 

gives  a  value  of  K  of 

K   =   0.014  Btu/hrft°F 
g 

and  this  is  in  close  agreement  with  published  experimental 
data  [34].   However,  if  Equation  12  is  evaluated  at  higher 
temperatures,  around  600°F,  its  theoretical  values  are 
about  18%  below  the  above  cited  experimental  value. 
This  discrepancy  arises  because  of  the  departure  of  real 
gases  from  the  assumptions  used  in  deriving  Equation  6. 

Table  6  has  been  prepared  to  show  a  relative  comparison 
for  air  and  argon  when  used  under  non-continuum  conditions 
and  at  typical  mean  temperatures  encountered  in  this 
study.   Table  6  suggests  that  the  smaller  particle  size 
inhibits  gaseous  conductivity  and  that  at  higher  tempera- 
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tures,  argon  performs  significantly  better  than  air, 
especially  in  the  pressure  range  of  interest  (0.50  ym  to 
10  mm  Hg) .   Carbon  dioxide  has  a  lower  thermal  conductivity 
than  air,  but  its  values  are  not  as  low  as  those  of  argon. 
In  checking  these  values  with  Equation  12,  keep  in  mind 
that  the  theoretical  conductivites  at  ambient  pressures 
tend  to  fall  as  much  as  20%  below  experimental  values,  as 
the  temperature  is  raised  from  32  to  600°F.   It  is  not 
apparent  if  this  deviation  continues  at  reduced  pressures. 

This  nominal  diameters  used  in  evaluating  the  powder 
particles  were  taken  as  the  average  of  the  screen  openings 
below  and  above  the  trapped  powder.   For  example,  if  powders 
were  trapped  between  Tyler  Sieve  #500,  screen  opening 
25  ym,  but  passed  through  Tyler  Sieve  #325,  screen  opening 
45  ym,  then  the  nominal  particle  size  was  taken  as  the 
average  of  the  two  values,  or  35  ym.   It  is  well  to  point 
out  that  the  nominal  diameter  is  used  for  convenience  in 
discussing  the  particles. 

Average  void  diameters  of  50%  of  the  nominal  diameter 
were  assumed  in  view  of  Figures  8  and  9  and  from  observation 
of  diatomaceous  earth  particles  under  an  optical  microscope. 
At  the  onset  of  this  study  it  was  assumed  that  powders 
separated  by  sieves  would  be  of  fairly  uniform  size. 
However,  Figures  8  and  9  indicate  that  the  nominal  size,  as 
defined  above,  is  larger  than  many  particles  found  in  the 
powders . 
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Particle- to -particle  conductivity.   Table  4  (page  80) 
shows  a  variety  of  vacuum  levels,  hot  face  temperatures, 
and  apparent  thermal  conductivity  values  for  powder  insula- 
tion under  vacuum.   These  values  were  obtained  from  typi- 
cal operating  conditions  of  a  high  temperature  thermal 
energy  storage  system.   With  the  experimental  data  of 
Table  5,  and  by  recalling  that  energy  transport  through 
powders  under  vacuum  occurs  principally  by  conduction,  one 
may  use  the  developed  theoretical  equations  to  proportion 
heat  flow  through  gas  and  particle-to-particle  conduction. 
This  will  be  done  in  the  following  paragraphs  by  evaluating 

the  conductivity  coefficients,  X     and  T   . 

g      s 

Table  7  illustrates  the  proportioning  of  apparent 
conductivity  into  gaseous  and  particle-to-particle 
conductivity.   The  gaseous  conductivity,  K  ,  was  calculated 
from  Equation  12  and  the  apparent  conductivity  values,  K  , 
were  taken  from  Table  5.   The  particle-to-particle 
conductivity  was  then  calculated  from  these  values  using 
the  estimated  area  ratios.   As  a  result  of  viewing  the 
diatomaceous  earth  powders  under  an  optical  microscope, 
a  parallel  heat  flow  model,  in  the  low  pressure  range,  was 
selected.   The  gas  and  solid  materials  were  viewed  as 
parallel  tubes.   Heat  flows  through  the  separate  tubes 
composed  of  different  materials,  in  parallel.   This  can 
be  expressed  as 
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TABLE  7 

Proportioning  of  Thermal  Conductivity 
Into  Gaseous  and  Solid-to-Solid  Conductivity 

Small  Grain  Particles 


Line* 

Tb 
(°F) 

P 
Cwm  Hg) 

K  * 
a 

(Btu/hrft°F) 

K 
g 
(Btu/hr£t°F) 

(Btu/hrft°F) 

1 

247.5 

500 

0.00364 

0.00204 

0.0316 

2 

417.0 

600 

0.00648 

0.00239 

0.0319 

3 

396.5 

800 

0.00750 

0.00286 

0.0364 

4 

209.5 

850 

0.00483 

0.00327 

0.0145 

5 

417.0 

1000 

0.00722 

0.00343 

0.0308 

6 

372.0 

1000 

0.00790 

0.00348 

0.0357 

7 

327.5 

1000 

0.00444 

0.00355 

0.0103 

8 

269.5 

1000 

0.00407 

0.00363 

0.0068 

9 

263.5 

1000 

0.00417 

0.00364 

0.0075 

10 

411.5 

760,000t 

0.0648 

0.01988 

0.3445 

11 

293.0 

760,000 

0.0460 

0.01850 

0.2172 

Large  Grain  Particles 

24 

268.5 

175 

0.00635 

0.00202 

0.0378 

25 

393.5 

200 

0.00656 

0.00211 

0.0388 

26 

348.0 

300 

0.00590 

0.00307 

0.0264 

27 

378.5 

800 

0.00723 

0.00642 

0.0131 

28 

212.0 

1000 

0.00972 

0.00758 

0.0252 

29 

308.0 

1000 

0.01198 

0.00750 

0.0445 

30 

404.0 

3800 

0.01343 

0.00739 

0.0572 

31 

326.0 

600 

0.01714 

0.01352 

0.0434 

From  Table  5. 

760,000  ym  Hg  equals  one  standard  atmosphere 


qT  =  V  +  %  <48) 


This  results  in, 

A 


K   =  -^  K   +  ^  K  (49) 

a    A~,   g   AT   s  v  J 

The  area  fractions  of  the  above  equation  would  be  a 
difficult  parameter  to  measure  for  a  loose  powder.   These 
values  can  be  estimated  under  the  assumption  that  the  area 
ratios  can  be  deduced  from  the  void  ratio.   For  example, 
if  all  of  the  loose  powder  material  in  a  unit  cube  were 
separated  into  a  uniform  solid  in  the  bottom  of  the  cube, 
with  the  void  volume  (residual  gas  space)  above  it, 
then  the  ratio  of  the  void  volume  to  total  volume  equals 
the  void  ratio.   Since  the  cube  is  of  unit  width,  the  void 
will  equal  the  area  ratio  as  given  in  the  first  term  on 
the  right-hand  side  of  Equation  49.   The  same  reasoning 
will  also  apply  to  the  solid  ratio  which  is  the  ratio  of 
the  solid  volume  to  total  volume  or  the  solid  area  to  the 
total  area,  as  seen  from  the  edge  of  the  unit  cube.   Hence, 
for  a  differential  thickness  of  powder,  the  projected 
solid  and  void  areas  will,  on  the  average,  be  equal  to  the 
solid  ratio  and  the  void  ratio,  respectively. 

In  reviewing  Table  5,  one  would  generally  expect  a 
relationship  between  particle-to-particle  conductivity  and 
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temperature.   This  relationship  certainly  exists,  but  is 
not  readily  apparent  as  it  is  masked  by  the  pressure 
sensitive  gaseous  conductivity.   This  is  evident  in  Lines 
4-5  of  Table  5.   These  values  suggest  the  need  for  more 
detailed  information  --  particle  sizes,  void  dimensions, 
and  so  forth.   Also,  the  temperature  gradient,  particularly 
near  the  hot  face  boundary,  influences  the  overall  apparent 
conductivity.   [Typical  temperature  profiles  are  shown 
in  Figure  22.]   For  amorphous  powder  particles  with  an 
increasing  value  of  particle  conductivity,  k  ,  the  conduc- 
tance next  to  the  hot  face  boundary  will  be  signifiantly 
more  than  throughout  the  remainder  of  the  material. 

One  particular  difficulty  in  working  with  powders 
is  that  the  value  of  its  solid  particle  conductivity,  k  , 
is  not  available.   An  obvious  choice  as  a  substitute  for 
k  of  diatomaceous  earth  would  be  to  use  the  conductivity 
coefficient  for  another  amorphous  material,  such  as  that 
of  pyrex  or  quartz  glass.   The  reason  for  selecting  these 
materials  is  that  the  electron  microscope  photographs 
reveal  particles  that  show  cleavage  planes  characteristic 
of  an  amorphous  material.   Evidence  of  crystallinity  was 
not  apparent.   Also,  the  data  reported  herein  indicate 
that  point-to-point  solid  conductivity  increases  with 
temperature,  and  this  is  consistent  with  amorphous  material 
such  as  pyrex  or  quartz  glass. 
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Equation  37  can  now  be  evaluated  for  particle-to- 
particle  conductivity.   Recall  Equation  37, 

Cks  ni  Dc  NA 
Ks  =   2Nh  ln(Dc/y*)  (37) 

where, 

,     (I  -  Mf)El  +  (1  .  y2)/E2   1/3 

?*    =    %  F[    1/Dl  +  1/D2 3)  (38) 

The  above  equations  will  be  evaluated  based  on  the  follow- 
ing assumed  values  for  a  one  foot  cube  sample  of  Tyler 
Sieve  #500  particle  (nominal  diameter  D   =  1.5  x  10"5  ft) 
at  a  mean  temperature  of  400°F. 

ks  =  0.76  Btu/hrft°F  (pyrex  glass) 

n   =  1  contact  point  per  particle  (heat  flow  in) 

Nh  =  1/I  Dn  "  -1--33  x  105  Particles/ft 

NA  =  (1/i  D2)2  -  1.78  x  1010  particles/ft2 

F   =  (2/3)(20.0  lbf/ft2)/1.78  x  1010  particles/ft2 

=  7.5  x  10" 10  lbf 
y  1  =  y2  =  0.25 
E2  =  E2  =  10  x  106  lb£/in2 

Dc  =  Dx  =  D2  =  3. 3  x  10"6  ft  (estimated  equivalent 
diameter  of  contact  region.) 

Note  in  the  above  values  that  the  particle  size  has 
been  estimated  as  1/2  the  nominal  particle  size.   The 
weight  of  force,  F,  is  taken  at  2/3  the  depth  of  the  powder 
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bed,  and  the  force  is  divided  by  the  number  of  points  of 
contact  that  support  it.   In  this  case  a  value  of  unity 
has  been  assumed.   This  will  be  explained  later.   The 
value  of  D  was  estimated  from  Figure  14  and  for  this 
example  all  particle  contact  points  are  assumed  to  have 
the  same  value. 

With  the  above  values,  Equation  37  gives 

Ks  =  C(0.09)  Btu/hrft2oF 

The  correction  factor  required  is  on  the  order  of  unity. 
This  is  consistent  with  values  of  solid  conductivity  shown 
in  Table  6  for  particle-to-particle  conductivity  values 
obtained  at  low  pressure. 

The  average  stress  holding  two  particles  in  contact 
for  the  conditions  listed  above  is  found  with  the  use 
of  Equations  21  and  35,  to  be 

S  =  2.1  x  104  lbf/in2 

This  value  appears  to  be  more  than  sufficient  to  maintain 
two  particles  in  intimate,  solid  contact.   It  is  on  this 
basis  that  the  contact  area  resistance  was  replaced  with 
a  contact  region  conductive  resistance. 

The  two  values  at  atmospheric  pressure  suggest  that 
other  factors  affect  heat  transfer  at  that  pressure  level. 
One  avenue  that  may  be  explored  is  that  at  conditions  where 
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Figure  14 

Electron  Microscope  Photograph  of  Surface 
of  Diatomaceous  Earth  Particle,  Tyler  Sieve  #500 
fx  10,00( 


Figure  15 


Electron  Microscope  Photograph  of  Surface 
of  Diatomaceous  Earth  Particle,  Tyler  Sieve  #170 
(x  10,000) 
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continuum  assumptions  apply,  Kn  <<  1,  a  series-parallel 
model  may  be  more  applicable.   Since  the  analysis  of  energy 
transport  through  powders  at  atmospheric  pressure  was 
beyond  the  scope  of  this  study,  this  topic  was  not  pursued. 

In  the  sample  calculation,  the  number  of  contact 
regions,  n,  was  assumed  to  be  unity.   This  value  was  selected 
after  inspecting  the  particles  in  an  optical  microscope  which 
gives  a  three  as  opposed  to  a  two-dimensional  view  as  seen 
in  the  electron  microscope  photographs.   In  three  dimensions 
the  diatomaceous  earth  particles  look  much  like  an  assembly 
of  tinker  toys.   The  particles  adhere  in  groups  and  these 
groups  are  separated  with  thin  columns  of  particles.   These 
columns  appear  like  small  magnets  clinging  together.   Within 
the  groups,  it  appears  as  though  the  individual  particles 
are  repelling  each  other.   This  separation  of  particles 
reduces  the  number  of  contacts  between  particles.   For 
example,  small  pebbles  form  piles  in  which  each  higher 
pebble  is  supported  by  several  lower  pebbles  with  several 
points  of  contact.   Because  of  the  formation  of  groups, 
separated  by  strands,  diatomaceous  earth  particles  seem 
to  minimize  the  number  of  contact  regions  without  breaking 
apart.   In  the  following  paragraphs,  the  cause  of  this 
arrangement  is  discussed. 

The  force  of  attraction  between  particles  is  thought 
to  be  electrostatic  in  origin.   In  that  sense,  it  is 
similar  to  polished  metals  which  adhere  to  each  other. 
The  diatomaceous  earth  particles  do  not  appear  to  interlock 
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as  do  pebbles  (uncharged  particles) .   It  was  thought  that 
water  surface  tension  would  be  a  factor  in  particle 
adhesion.   However,  heating  the  particles,  to  drive  off 
moisture,  did  not  appear  to  alter  their  attraction,  at 
least  at  atmospheric  pressure. 

When  the  diatomaceous  earth  particles  were  viewed 
under  an  optical  microscope  and  the  glass  support  was 
vibrated  by  tapping,  the  individual  particles  did  not 
vibrate.   This  is  in  contrast  to  similar  observations  of 
aluminum  powder  particles,  which  have  no  electrostatic 
forces.   The  absence  of  induced  vibrations  implies  a  force 
of  attraction  between  the  particles. 

One  can  see  from  Equation  35  that  the  effect  of  a  force 
or  weight  on  porous  powders  is  two  fold.   First,  it  increases 
the  contact  diameter.   This  effect  is  small  as  evidenced  by 
the  log  term  in  the  denominator.   Second,  and  more  signifi- 
cantly, increased  force  presses  more  particles  together 
and  as  a  result  a  larger  number  of  contact  regions  form. 
This  net  effect  of  weight  on  K  ,  as  shown  in  Figure  16, 
contrasts  the  curves  for  Tyler  Sieve  #500  powders  for  com- 
pressed and  uncompressed  material.   Note  that  the  effect  of 
compressing  the  powder  at  lower  temperatures  is  to  decrease 
its  thermal  conductivity  as  compared  to  uncompressed  particles 
Presumably,  this  occurs  because  of  a  decrease  in  the  thermal 
conductivity  of  the  gas  that  results  from  a  reduced  void 
dimension  in  the  non-continuum  regime. 
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At  higher  temperatures,  on  the  other  hand,  the  apparent 
conductivity  of  the  compressed  powders  increases  rather 
sharply.   Nevertheless,  Figure  16  suggests  that  compressed 
powders  may  be  used  to  support  a  storage  tank  weight  without 
a  large  penalty  in  insulation  performance.   Specifically, 
fine  compressed  powders  can  be  placed  under  a  container  for 
support  and  loose,  uncompressed  powders  can  be  placed 
in  the  remaining  volume.   In  this  manner  only  limited  supports 
or  no  supports  at  all  are  needed  to  separate  the  internal 
and  external  tank.   This  eliminates  the  need  for  large 
internal  or  external  supports  through  which  heat  may  be 
lost. 

The  apparent  thermal  conductivity  values  of  Table  5 
have  been  reviewed  except  for  the  last  four  categories: 
raw  diatomaceous  earth,  pulverized  diatomaceous  earth,  air 
only,  and  fiberglass  cloth  and  aluminum  foil.   The  raw 
diatomaceous  earth  appears  to  perform  similar  to  the  powder 
in  the  #170  pan.   That  is,  the  powder  particles  and  the 
average  void  dimensions  are  large.   Pulverizing  the  raw 
diatomaceous  earth  and  using  it  without  separating  it  into 
different  particle  sizes  will  be  of  some  advantage. 

The  air-only  condition  was  measured  to  demonstrate 
the  effect  of  using  air  alone  without  convection  suppression 
or  radiation  shielding. 

The  fiberglass  cloth  and  aluminum  foil  demonstrate  a 
value  that  one  might  expect  from  this  material.   The  one 
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value  given  in  Table  4  is  relatively  high,  but  it  can  be  made 
some-what  lower  by  using  more  layers  of  fiberglass  cloth 
and  aluminum  foil.   However,  note  that  this  multilayer 
material  is  at  a  pressure  of  0.40  mm  Hg.   To  produce  a 
significant  drop  in  K   for  this  material,  the  pressure  must 
be  reduced  to  a  level  of  0.10  mm  Hg  or  less.   See  Table  5, 
page  79 . 

The  fiberglass  cloth  and  aluminum  foil  were  difficult 
to  put  together  since  the  multi-foil  layers  were  individually 
wrapped.   Some  of  the  probles  with  wrapping  the  insulation 
were 


1.  the  need  for  securing  the  layers  so  that  they 
did  not  unravel, 

2.  a  support  system  was  required  to  separate  the 
inner  and  outer  tank,  and 

3.  the  supports  and  piping  to  and  from  the  inner 
tank  had  to  be  carefully  wrapped  to  prevent 
heat  loss. 


It  is  well  to  point  out  that  multi-foil  insulation  at  low 
pressure  works  well.   But,  its  complexity  and  expense  make 
it  an  unlikely  candidate  for  routine  high  temperature  use. 

Effect  of  particle  size  in  a  non-continuum.   Figure 
16  (page  96)  illustrates  the  trend  of  particle  size  versus 
conductivity.   The  larger  particles,  in  the  case  from  Tyler 
Sieve  #170,  pass  heat  about  three  times  more  readily  than 
the  particles  of  Tyler  Sieve  #500.   This  trend  is  also  seen 
in  the  cryogenic  region  [10].   The  apparent  reason  for  this 


is  that  small  powders  suppress  gaseous  conductivity  more 
effectively  and  they  create  a  greater  resistance  to  heat 
flow  due  to  contact  region  resistance  as  explained  in 
Chapter  3.   It  is  interesting  to  note  the  behavior  of  the 
compressed  powder  (Tyler  Sieve  #500)  in  Figure  16.   At 
the  low  temperature  end  of  the  graph  the  compressed  powders 
have  a  high  resistance  to  heat  flow  due  to  non-continuum 
conditions  created  by  the  fine  particles.   However,  at 
higher  temperature,  the  conductivity  increases  rapidly.   It 
is  assumed  that  this  is  due  to  the  presence  of  additional 
solid  material  through  which  heat  may  pass.   This  powder 
was  compressed  to  about  one  half  of  its  original  volume. 
The  above  observation  is  consistent  with  Equation  37  in  which 
heat  flow  is  directly  proportional  to  the  number  of  particles, 
NA. 

Effect  of  vacuum  level.   The  vacuum  level  has  a  major 
impact  on  the  heat  loss  characteristic  of  the  powder  insulation 
under  vacuum.   Figure  17  is  a  theoretical  curve  showing  the 
characteristic  decline  of  gaseous  conductivity  for  air  and 
argon.   The  decline  is  seen  to  be  linear  at  points  less  than 
1000  ym  Hg.   The  gaseous  conductivity  below  this  pressure 
varies  asymptotically  with  pressure.   Note  the  relatively 
flat  response  at  pressures  near  10,000  ym  Hg  (10  mm  Hg)  in 
Figure  17.   This  indicates  the  importance  of  operating 
evacuated  insulation  at  less  than  10  mm  Hg. 

Figure  17  illustrates  the  trend  in  the  gaseous  con- 
ductivity for  air  and  argon  at  pressures  below  10  mm  Hg 
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where  a  rapid  declind  in  conductivity  occurs.   On  this 
scale,  a  recommended  operating  range  falls  between  0.5 
mm  Hg  to  3  mm  Hg  since  these  values  are  easily  obtained 
with  common  mechanical  vacuum  pumps. 

Selection  of  a  particular  operating  pressure  will 
depend  upon  the  particle  size,  particle  material,  residual 
gas  present,  temperature  level,  cost  of  stored  energy  and 
so  forth.   A  pressure  much  lower  than  0.5  mm  Hg  can  be  used 
but  at  the  expense  of  more  complex  vacuum  pumps  (diffusion 
or  molecular  pumps),  longer  pumping  times,  and  potential 
problems  of  outgassing.   From  theoretical  considerations  a 
point  will  be  reached  in  which  further  reductions  in  heat 
loss  through  gas  conduction  becomes  negligible  compared 
to  solid  particle  conduction  losses.   This  investigation 
did  not  experimentally  identify  this  pressure  region  due  to 
equipment  limitations.   However,  part  of  the  goal  of  this 
work  was  to  identify  a  practical  upper  limit  for  pressure 
commensurate  with  a  low  coefficient  of  thermal  conductivity 
(K  <  0.015  Btu/hrft°F  (@  1000°F).   The  pressure  range  of 
0.5  mm  Hg  to  3.0  mm  Hg  is  designated  as  an  upper  pressure 
range  for  accomplishing  this  goal. 

Effect  of  carbon  and  iron  in  diatomaceous  earth.   Carbon 
was  added  to  the  diatomaceous  earth,  at  about  1%   by  weight, 
to  enhance  pouring  and  separating  characteristics  of  the 
latter.   Very  small  amounts  of  steel  powder  found  its  way 
into  the  material  that  was  pulverized  in  a  ball  mill.   [The 
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ball  mill  container  and  rollers  were  made  of  steel.]   From 
the  limited  data  of  Table  5,  it  appears  that  neither 
the  carbon  nor  fine  steel  had  a  significant  measurable 
effect  on  the  thermal  conductivity  of  diatomaceous  earth. 

Thermal  Energy  Storage  System 

The  thermal  energy  storage  system  was  filled  with 
approximately  500  pounds  of  diatomaceous  earth  that  had 
been  pulverized  in  a  ball  mill  for  30  minutes.   Approxi- 
mately 1%   by  weight  of  carbon  black  was  added  prior  to 
pulverizing  since  the  carbon  was  found  to  give  the  dia- 
tomaceous earth  better  pouring  characteristics.   A  small 
amount  of  raw  diatomaceous  earth,  50  pounds,  had  to  be  used 
to  completely  fill  the  tank.   This  material  was  located  in 
the  top  back  section  of  the  tank  and  no  observable  tempera- 
ture difference  in  this  area  was  noticed  during  the  test. 

The  tank  was  filled  by  first  evacuating  the  insula- 
tion space  and  then  letting  the  powders  be  pushed  into 
the  tank  as  atmospheric  pressure  filled  the  annular  void. 
This  procedure  was  much  quicker  and  easier  than  trying  to 
fill  the  insulation  space  by  hand.   It  was  also  found  that 
better  packing  and  distribution  of  the  insulation  material 
was  obtained  since  no  hot  spots  were  noticeable  on  the 
exterior  of  the  outer  tank. 

The  storage  tank  had  thermocouples  attached  at  various 
points  to  the  exterior  of  the  inner  tank,  as  shown  in 
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Figure  18,  and  they  were  also  placed  in  the  powder  insula- 
tion so  that  a  temperature  profile  through  the  powder  could 
be  observed.   Surface  temperatures  of  the  outer  tank  were 
taken  externally. 

The  thermocouple  port  was  part  of  the  forward  vacuum 
port.   The  thermocouples  were  fed  through  a  u-shaped 
tube  and  the  u-shaped  tube  was  then  filled  with  a  high 
temperature  wax.   This  formed  an  effective  seal  since  the 
thermocouple  port  was  far  enough  away  from  the  heated  tank 
so  that  the  thermocouples  were  not  affected  by  heat  conduc- 
tion through  the  wires.   That  is,  heat  did  not  cause  the 
wires  to  expand  and  contract.   As  a  result,  no  leaks  were 
detected  between  the  wire  and  wax  seal.   The  thermocouples 
were  stripped  of  insulation  prior  to  emerging  from  the 
wax  to  prevent  leaks  from  developing  under  the  wire  insula- 
tion. 

The  cylindrical  portion  of  the  inside  of  the  outer 
tank  was  lined  with  a  blanket  of  fine  steel  wool,  grade 
(000)  .   This  provided  the  path  of  least  flow  resistance  for 
evacuating  gases  from  the  powder  insulation.   It  appears  that 
the  steel  wool  made  it  possible  for  air  to  flow  essentially 
in  a  radial  path  out  of  the  powders  so  that  an  air  molecule 
traveled  approximately  6.75  inches  through  the  powder 
instead  of  through  all  of  the  powder  from  the  bottom  to  top 
of  the  tank.   It  also  filtered  powder  entrained  in  the  evacu- 
ation flow  so  that  only  minimum  filtration  was  needed  at  the 
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inlet  to  the  vacuum  pump.   The  steel  wool  worked  well  because 
it  slowed  the  flow  of  air  sufficiently  so  that  entrained 
particles  could  not  be  carried  any  great  distance.   This 
was  observed  in  a  glass  cylinder  that  was  closed  at  its 
bottom  and  evacuated  at  its  top. 

The  steel  wool  blanket  was  fabricated  by  weaving  the 
4"  wide  wool  strands  through  wire  fencing.   The  fencing  was 
rolled  up  and  inserted  into  the  outer  tank  where  the  fencing 
expanded,  holding  the  steel  wool  against  the  tank  surface. 

The  inner  tank  was  charged  with  180  pounds  of  vegetable 
oil  (soybean) .   The  powders  were  evacuated  and  the  tank 
heater  was  turned  on  and  heating  of  the  oil  was  initiated. 
At  intervals  the  tank  was  allowed  to  reach  steady-state,  so 
that  the  effectiveness  of  the  insulation  could  be  determined. 
When  the  oil  reached  a  temperature  near  its  boiling  point, 
the  power  was  disconnected  and  the  tank  was  allowed  to  cool. 

Results  of  Prototype  Thermal  Energy  Storage  System  Test 

A  typical  result  of  the  cool  down  process  is  shown  in 

Figure  19  for  a  test  with  the  powder  evacuated  and  a  test 

in  which  the  vacuum  on  the  powder  was  released  after  the 

tank  of  oil  was  heated  close  to  its  boiling  point  (T,  -  -,  . 

&  r      v  boiling 

=  700°F  @  1  atm.  for  new  oil)  [3]. 

Figure  20  shows  typical  temperature  profiles  through  the 
evacuated  powder  insulation  for  several  insulation  hot  face 
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600- 


oA  -  curve  for  8/12/80 
11:30  am 

*B  -  curve  for  8/16/80 
8:00  am 

C  -  curve  for  8/19/80 
9:00  am 


Average  ambient  tempera- 
ture 84°F 
(Data  from  Appendix  A) 


Distance  From  Inner  Tank 
Through  Insulation  (inches) 


Figure  20 

Temperature  Profile  Through 
Powder  Insulation  Under  Vacuum 
(2-3  mm  Hg) 
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temperatures.   The  data  for  Figures  19  and  20  were  taken 
from  the  data  presented  in  Appendix  A. 

All  thermocouples  were  calibrated  to  read  within 
1/2°F  of  each  other  when  immersed  at  212°F  and  at  70°F. 
When  all  thermocouples  were  fixed  in  place,  they  were 
checked  again  at  room  temperature  to  insure  uniformity. 

An  electric  resistance  heater  was  built  into  the 
storage  tank  as  shown  in  Figure  7(a).   It  consisted  of 
a  glass  rod  with  a  resistance  heating  coil  wound  around 
it  and  the  wire  was  covered  with  a  thin  asbestos  cover, 
for  electrical  insulation.   The  heater  was  placed  in  the 
heater  port.   A  rheostat  was  used  to  regulate  AC  power 
to  the  heating  element.   An  ammeter  and  a  voltmeter  were 
connected  to  the  heating  element  leads  to  measure  the 
power  delivered  to  the  tank.   See  Figure  21. 

Discussion  of  Results  for  the  Thermal  Energy  Storage  System 
Test  ~"  ' ' 

Insulation  performance.   Figure  19  indicates  that 

powder  insulation  under  vacuum  performed  as  expected  as  a 

high  temperature  thermal  insulation.   Using  70 °F  as  the 

reference  for  the  sensible  energy  obtainable  from  the 

storage  tank,  Figure  19  shows  a  4.0%  loss  of  the  stored 

energy  the  first  day.   After  one  week,  over  50%  of  the 

original  stored  energy  was  still  available. 
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In  contrast  to  the  evacuated  insulation,  atmospheric 
insulation  (curve  II  in  Figure  19)  lost  approximately  60% 
of  its  energy  in  less  than  one  day,  once  the  vacuum  was 
released.   Specifically,  in  less  than  18  hours,  its  tempera- 
ture dropped  to  a  level  that  took  the  evacuated  powder 
insulation  168  hours  (7  days)  to  reach. 


Temperature 
Recorder 


Vacuum 
Gauge 


Vacuum 
Pump 


Storage  System 


amp 


z 


o 


-c 


Rheostat 


volts 


Figure  21 
Storage  System  Test  Equipment  Arrangement 


Figure  20  shows  typical  temperature  profiles  through 
the  powder  insulation  under  vacuum.   The  data  points  indica- 
ted on  the  graph  were  taken  by  placing  thermocouples  in 
the  powder  at  1 ,  3  and  5  inches  from  the  hot  surface  boundary 
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The  cold  face  temperature  was  measured  on  the  surface  of 
the  external  tank.   The  thermocouple  locations  of  the  data 
points  of  Figure  20  are  shown  in  Figure  13  and  a  reference 
list  is  shown  in  Appendix  A. 

The  temperature  profiles  of  Figure  20  raise  some  inter- 
esting questions  concerning  the  reasons  for  the  rapid 
increase  in  the  slope  of  the  curve  as  the  hot  surface 
temperature  increases  and  as  the  hot  surface  boundary  is 
approached.   Several  causes  for  this  are 

1.  for  cylinderical  geometry  where  energy  transport 
is  radial,  a  decreasing  temperature  profile  is 
expected  since  the  amount  of  insulating  material 
increases  with  an  increasing  radius, 

2.  the  particle  thermal  conductivity  may  decrease 
with  a  rise  in  temperature, 

3.  the  powder  next  to  the  hot  surface  boundary  may 
impede  radiative  transfer, 

4.  the  local  effect  of  temperature  on  lowering 
gaseous  molecular  energy  transfer  is  more 
pronounced  near  the  hot  surface  boundary  than 
near  the  cold  surface  boundary,  and 

5.  an  interaction  between  the  above  mechanisms  may 
be  occurring  with  a  consequential  local  rise  in 
temperature  as  the  hot  surface  is  approached. 

The  shape  of  Curve  A  in  Figure  20  suggest  that  radiant 
and  gaseous  molecular  energy  transfer  may  be  responsible 
for  the  temperature  gradient  next  to  the  hot  surface 
boundary.   This  is  evidenced  by  the  sudden  rise  in  tempera- 
ture of  Curve  A  between  a  point  three  inches  from  the  hot 
surface  and  the  hot  surface.   Hence,  although  the  net 
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effect  of  radiant  transport  through  the  powders  may  be 
negligible  its  local  effect,  in  combination  with  the  other 
modes  of  energy  transport  near  the  hot  surface  boundary 
is  significant.   The  shape  of  the  right-hand  portion  of 
Curve  A,  as  well  as  Curve  B  and  Curve  C,  suggest  that 
linear  heat  transfer  mechanisms  control  energy  transfer 
in  that  region.   Also  that  portion  of  the  curves  are  more 
consistent  with  constant  heat  flux  distribution  through 
cylindrical  geometry  according  to  the  term  In  (r  /r.)  in 
Equation  41. 

The  temperature  profiles  shown  in  Figure  20  indicate 
that  powder  insulation  under  vacuum  is  effective  since 
the  outer  boundary  temperature  was  only  a  few  degrees 
higher  than  ambient. 

Heating  of  organic  oil.   Curve  I  of  Figure  19,  page  106, 
begins  at  a  temperature  well  below  the  normal  boiling  point 
of  soybean  oil,  that  is  700°F  at  atmospheric  pressure.   In 
the  first  run  of  the  high  temperature  storage  system,  the 
oil  was  heated  to  the  vicinity  of  its  boiling  point.   The 
consequent  rise  in  tank  pressure  caused  tank  failure. 
This  was  not  anticipated  as  the  storage  system  had  been 
pressure  tested  to  250  psig  and  the  air  and  oil  vapor 
pressure,  at  an  oil  temperature  of  700°F,  was  expected  to 
be  less  than  75  psig.   Analysis  of  the  tank  failure  showed 
that  the  heater  area  was  undersized.   Since  the  expected 
natural  convection  coefficient  of  heat  transfer  of  organic 
oil  is  small,  about  7.5  Btu/hrft2oF  at  670°F,  the  tempera- 
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ture  of  the  oil  next  to  the  heater  surface  can  reach  or 
exceed  its  boiling  point.   This  indicates  that  if  organic 
oil  is  used  for  heat  transfer  purposes,  then  the  heat 
transfer  surface  must  be  large  or  the  temperature  differences 
between  it  and  the  oil  must  be  small,  if  the  oil  is  to  be 
used  near  its  boiling  temperature.   This  also  implies  that 
for  forced  convection  systems,  high  flow  rates  will  be 
required  [3].   Also,  if  the  oil  is  to  be  used  near  its 
boiling  point,  then  safety  precautions  must  be  taken  in 
view  of  the  high  pressures  that  may  be  encountered. 

Steel  wool  blanket.   The  steel  wool  blanket  (grade  000) 
performed  well  in  suppressing  powder  entrainment  during 
the  evacuation  of  air  from  the  powder  insulating  material. 
The  operating  principle  of  the  steel  wool  blanket  was 
observed  by  evacuating  powder  in  a  glass  cylinder.   As 
the  air  was  evacuated,  small  holes  formed  in  the  powder 
permitting  air  to  escape.   As  the  air  exists  through  the 
small  holes,  its  velocity  in  the  cavity  was  sufficient  to 
entrain  powder  particles.   When  steel  wool  was  placed 
over  the  powder  surface,  the  air  jets  coming  out  of  the 
powder  surface  were  diffused  by  the  steel  wool,  thus  mini- 
mizing entrainment. 

Overall  storage  vessel  performance.   During  the 
heating  phase  of  the  storage  vessel  performance  test  the 
system  was  periodically  allowed  to  attain  a  steady  state 
condition.   At  that  point  the  current  and  voltage  to  the 
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heater  element  were  monitored  and  the  apparent  conductivity 
of  the  insulation  material  was  approximated  by  the  following, 


«loss  "  Ka  *  4T  Cln^o/Lr,)+  rat)  (50) 


The  first  part  of  Equation  50  represents  the  heat  loss 
form  the  cylindrical  portion  of  the  tank  and  the  second 
part  represents  the  heat  loss  through  the  ends  of  the  tank, 
The  diameter  of  the  end  of  the  tank  was  approximated 
by  taking  the  arithmetic  mean  of  the  inner  and  outer  tank 
diameters  to  compensate  for  edge  losses.   Solving  Equation 
50  for  the  apparent  conductivity,  K  ,  and  evaluating  with 

a. 

the  following  conditions, 

Q.   c  =  200  Btu/hr 
^■loss  ' 

rQ  =  1.25  ft  (outer  tank  radius) 

r^  =  0.67  ft  (inner  tank  radius) 

L  =  4.17  ft  (inner  tank  length) 

AT  =  (530  -  93)°F 

AX  =  0.58  ft  (insulation  thickness) 

D  =  1.92  ft 
gives, 

K   -  0.0097  Btu/hrft°F 

This  value  is  consistent  with  apparent  conductivity  values 
obtained  in  the  conductivity  test  cell.   It  should  be 
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recalled  that  the  density  of  the  powder  in  the  storage 
vessel  was  about  32  lbm/ft3.   With  this  density,  the 
vacuum  level  was  allowed  to  operate  between  2-3  mm  Hg. 

Heat  losses  through  the  inlet  and  outlet  pipes,  the 
thermocouple  wires,  and  the  tank  support  cables  were  found 
to  be  negligible  as  explained  in  Chapter  V.   The  heater 
port  was  found  to  carry  about  5%    of  the  total  heat  loss. 
If  a  larger  heater  port  should  be  needed,  as  indicated 
earlier,  then  heat  loss  through  it  could  become  significant, 
and  care  should  be  taken  in  its  design. 

Tank  support.   The  small  stainless  steel  cable  used 
in  the  high  temperature  thermal  energy  storage  vessel 
was  found  to  be  sufficient  for  its  intended  purpose.   Heat 
loss  through  the  calbes  was  found  to  be  negligible  as 
indicated  in  Chapter  V. 

The  test  data  for  the  Tyler  Sieve  #500  compressed 
powder,  Figure  16,  indicates  that  it  may  be  possible  to 
let  the  insulating  powder  support  the  inner  tank.   This  could 
be  accomplished  by  using  compressed  powders  under  the  tank 
and  then  filling  the  remaining  volume  with  loose  fill  or 
uncompressed  powders.   The  benefit  of  this  design  would  be 
the  elimination  of  a  heat  loss  path  and  the  elimination  of 
external  supports  needed  to  support  the  weight  of  the  inner 
tank  and  its  heat  storage  medium. 

Cost  of  high  temperature  storage.   Inevitably,  cost 
will  be  a  primary  factor  in  governing  the  use  of  the  high 
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temperature  thermal  energy  storage  systems.   Assume  a  self- 
sufficient  residential  system  requires  a  storage  capacity 
of  106  Btu.   The  required  storage  volume  is  then  approxi- 
mated as  follows 


Heat  Storage 

Material  Volume  (ft3)    Sphere  Radius  (ft) 

Soybean  Oil  (600°F)        65  2.5 


The  oil  storage  tank  will  require  two  tanks.   The  outer 
tank  radius  is  assumed  to  be  0.5  feet  larger  than  the 
inner  tank  radius.   Table  8  gives  an  estimated  cost  of 
this  storage  system. 

The  economics  of  the  high  temperature  storage  vessel 
are  greatly  enhanced  if  the  latent  heat  of  fusion  of  a 
heat  storage  medium  is  used  along  with  sensible  heat 
storage.   For  example,  if  NaN03  [35]  were  used  instead  of 
organic  oil,  then  the  high  temperature  storage  vessel  would 
contain  approximately  2.13  x  106  Btu.   This  gives  a  storage 
cost  of  $493  -  1221  per  106  Btu.   This  phase  change  material 
requires  a  rather  elaborate  heat  exchanger  and  it  will 
displace  some  of  the  heat  storage  material.   This  will  make 
the  above  cost  higher  but  the  cost  is  competitive  with  that 
of  low  temperature  thermal  energy  storage  vessels  such  as 
hot  water  storage. 
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Table  8 
Estimated  Manufacture's  Cost  for  Storage  Vessels 

System  Specifications 

Tank  volume  (ft3) 

(a)  inner  65 

(b)  outer  113 

Insulation  volume  (ft3)  48 

Tank  material  -  carbon  steel 

Cost  Range  of  Fabrication  [36] 

Spherical  tanks 

(two  required)  $500  -  2000+ 

Insulation  powder  $  50  -   100* 

(diatomaceous  earth) 

Labor  -  40  hrs  $400 

§  $10.00/hr 

Tank  testing  -  5  hrs  $100 

@  $20/hr 


Cost  for  storing  106  Btu  $1050  -  2600 

(Storage  medium  -  soybean 
oil  @  600°F) 


Cost  depends  upon  quantity ,  source  of  purchase  and  wall 
thickness  of  tanks. 

* 

Cost  depends  upon  quantity,  grade  of  material,  particle 
size  and  source  of  purchase. 
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Uncertainty  in  the  Measurement  of  K 
[ a 

The  uncertainty  analysis  is  presented  in  Appendix  C 
concludes  that  a  maximum  uncertainty  in  the  physical 
measurements  of  the  apparent  thermal  conductivity  coefficient 
to  be  approximately  4.11  at  a  specified  pressure.   This 
value  is  within  acceptable  engineering  limits  for  experi- 
mental work.   Included  in  Appendix  C  is  also  a  discussion 
concerning  the  methodology  in  the  use  of  thermocouples  for 
temperature  measurement. 


CHAPTER  VI 
CONCLUSIONS  AND  RECOMMENDATIONS 


Conclusions 


Based  on  the  work  described  in  this  dissertation, 
the  following  conclusions  are  given. 


1.  The  overall  use  of  powder  insulation  under 
vacuum  for  high  temperature  thermal  energy 
storage  is  technically  and  economically 
feasible. 

2.  Thermal  conductivity  values  of  0.0029  -  0.110 
Btu/hrft°F  were  achieved  experimentally  for 
vacuum  levels  ranging  from  0.5  -  3.0  mm  Hg  and 
hot  surface  temperatures  ranging  from  340  - 
711°F. 

3.  A  vacuum  level  of  0.5  mm  Hg  to  10  mm  Hg  is  the 
recommended  pressure  range  for  operating  powder 
insulation.   The  choice  of  vacuum  level  will 
depend  upon  particle  size,  particle  material, 
type  of  gas  present,  temperature  level,  and 
system  costs. 

4.  Small  particles  of  diatomaceous  earth  suppress 
gaseous  conductivity  more  effectively  than  large 
particles  when  compared  using  the  same  residual 
gas  of  the  same  pressure. 

5.  From  the  theoretical  considerations  the  apparent 
thermal  conductivity  of  powder  under  vacuum  can 
be  improved  by  selecting  residual  gases  that 
inhibit  molecular  energy  transport,  such  as  gases 
with  high  molecular  weight  and  specific  heat 
ratios . 

6.  An  evacuated  powder  can  be  used  to  support  a  load 
and  still  maintain  a  relatively  low  thermal  con- 
ductivity. 
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The  desirable  characteristics  of  an  insulating 
powder  particle  are  small  particle  size,  spheri- 
cal geometry,  low  thermal  conductivity,  non-con- 
ducting electrically,  non- corroding,  stable  at 
elevated  temperatures,  universally  available,  and 
inexpensive. 

Particle  size  and  distribution,  as  measured  by 
mechanical  spearators,  of  the  materials  used  in 
this  study  did  not  appear  to  be  consistent  with 
electron  microscope  photographs  of  the  material. 
The  particles  as  viewed  by  the  electron  micro- 
scope appear  to  be  much  smaller  and  of  various 
sizes . 

The  use  of  a  fine  steel  wool  blanket  permitted 
evacuation  of  air  from  the  storage  vessel  with 
minimal  entrainment  of  the  powder  insulation. 


Recommendations 

Based  on  the  work  carried  out  in  developing  this 
dissertation,  the  following  recommendations  are  made. 

1.  There  is  need  for  quantitative  data  and  correla- 
tion techniques  for  predicting  actual  particle 
size,  shape,  and  distribution. 

2.  Methods  can  be  explored  for  fabricating  porous 
powders  that  will  minimize  heat  transfer  by 
particle-to-particle  conductivity  while  at  the 
same  time  suppressing  gaseous  and  radiative 
conductivity. 

3.  Experimental  data  at  low  pressure  and  high 
temperature  for  apparent  thermal  conductivity 
of  porous  powders,  with  various  interstitial 
gases,  would  be  very  useful. 

4.  The  evaluation  of  the  temperature-pressure 
relationship  of  organic  oils  along  with  an 
evaluation  of  the  corrosion  potential  of  the 
oils  at  elevated  temperatures  would  be  useful 
in  designing  high  temperature  thermal  energy 
storage  vessels. 
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There  is  a  need  for  an  investigation  into  the 
design  of  high  temperature  storage  vessels  for 
use  with  phase  change  material.   The  primary 
area  of  investigation  would  be  to  identify  and 
develop  suitable  heat  exchangers  or  other  means 
of  extracting  energy  from  the  solidifying  phase 
change  material. 


APPENDIX  A 
STORAGE  VESSEL  TEST  DATA 


Thermocouple  Number  and  Location 
For  Prototype  Thermal  Energy  Storage  Vessel  Test 


Thermocouple 

No.  * 

Location  (See 

Figure  18) 

1 

Bottom,  front 

2 

Bottom  middle 

3 

Bottom,  back 

Positions  on 

4 

Top,  back 

inner  tank 

5 

Top,  middle 

1  -  7 

6 

Top,  front 

7 

Center,  front 

8  Bottom  Inlet  Pipe  4-1/2"  from 
inner  tank 

9  Heater  Port  4-1/2"  from  inner  tank 

10  Top  Inlet  4-1/4"  from  inner  tank 

11  Bottom,  front  4"  from  inner  tank 

12  Bottom,  back  4"  from  inner  tank 

13  Side,  front  5"  from  inner  tank 

14  Side,  front  3"  from  inner  tank 

15  Side,  front  1"  from  inner  tank 

16  Side,  back  4"  from  inner  tank 

17  Top,  back  4"  from  inner  tank 

18  Top,  front  4"  from  inner  tank 

19  Center,  back  4"  from  inner  tank 

20  Heater  Port  -  on  outside  of  heater 
rod  in  heater  port 

21  Outer  tank  surface  (moveable  probe) 
2  2  Ambient 


Chromel-Alumel  thermocouples,  referenced  to  ice  bath. 
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APPENDIX  B 

CONVECTION  HEAT  TRANSFER  ANALYSIS 
FOR  A  POROUS  POWDER 


In  powder  with  particle  dimensions  less  than  0.2  inches 
in  diameter,  investigators,  such  as  Hill  and  Wilhelm  [20], 
found  no  evidence  of  convective  heat  transfer  in  quiescent 
gas-solid  beds.   The  condition  encountered  in  their  study- 
was  similar  to  those  of  this  work. 

The  above  investigation  applied  particularly  to  a 
gas-solid  bed  at  atmospheric  pressure.   This  study  was 
primarily   interested  in  heat  transfer  through  powders 
under  vacuum  conditions.   Hence,  for  the  conditions  mentioned 
above,  it  was  concluded  that  convective  heat  transfer 
would  be  eliminated.   Justification  for  this  conclusion 
follows . 

Convective  heat  transfer  in  the  laminar  region  can 
be  expressed  in  the  following  form, 


Nu  =  (j£)  =  G(Ra)%  (51) 

Kf 


where  b  is  a  characteristic  dimension  of  a  given  geometry, 
G  is  a  constant  of  the  order  of  unity  and  Ra  is  the  Rayleigh 
number.   In  expanded  form  Equation  51  can  be  written  as, 


Nu  =  G(gBCpp2b3AT/yk£)^  (52) 


The  term  g  for  a  perfect  gas  is  3  =  1/T.   The  task  now  is 
to  show  that  convective  heat  transfer  is  negligible  for 
conditions  encountered  in  this  study.   Specifically,  one 
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must  show  that  convective  heat  transfer  is  negligible 
over  the  temperature  and  pressure  range  encountered. 
These  values  are  250-1000°F  and  0.5-760  mm  Hg. 

If  Equation  52  is  evaluated  at  the  250°F  and  1000°F 
at  a  pressure  of  one  atmosphere,  for  b  =  0.015  ft,  AT  = 
10°F  ,  the  Nusselt  number  is  found  to  be  of  the  order  of 
unity  for  both  temperatures.   Since  the  Nusselt  number 
is  given  by  Equation  2,  then 


Nu  =  !*  -  1  (2) 


and  the  convective  heat  transfer  coefficient  is  seen  to  be 
of  the  same  order  of  magnitude  as  the  thermal  conductance 
of  the  fluid.   That  is, 

h  -  k/b  (53) 

A  decrease  in  the  heat  transfer  coefficient  with  a 
decrease  in  pressure,  as  calculated  by  Equation  52,  is 
illustrated  by  substituting  for  the  density  term, 

p  =  P/RT  (7) 


Assumed  maximum  temperature  difference  between  particles 
or  between  particles  and  a  container  boundary. 
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Nu  =  G(gC  P2b3AT/y£kfR2T2)%  (52) 

Restricting  the  above  equation  to  a  continuum,  it  is  clear 
that  the  Nusselt  number  and,  hence,  the  convective  heat 
transfer  coefficient  decreases  as  the  pressure  decreases. 
The  decrease  is  proportional  to  P2.   More  importantly,  the 
above  equation  shows  that  the  Nusselt  number  is  proportional 
to  (P/T)  2.   Typically,  the  operating  conditions  encountered 
in  this  study  were  that  pressure,  P,  would  be  decreased  and 
temperature,  T,  would  increase.   Assuming  a  continuum 
exists  in  a  gas-solid  bed  for  a  fluid  at  1  mm  Hg,  and  temera- 
ture  increased  from  250°F  (710°R)  to  1000°F  (1460°R),  the 
Nusselt  number,  and  hence  the  convective  heat  transfer 
coefficient,  would  decrease  by  approximately  two  orders 
of  magnitude,  as  compared  to  values  obtained  at  atmospheric 
conditions.   The  above  factors  indicate  that  convective 
heat  transfer  is  no  longer  significant  at  reduced  pressures 
and  high  temperatures. 

The  above  paragraph  illustrates  that  convective  heat 
transfer  becomes  negligible  as  the  pressure  changes  from 
760  mm  Hg  to  1  mm  Hg.   In  the  region  of  1  mm  Hg  and  where 
the  gas  molecules  are  confined  such  that  Kn  >  1 ,  the  assump- 
tion of  a  continuum  no  longer  applies.   In  this  case,  convec- 
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tion  heat  transfer  is  negligible  and  molecular  thermal 
conductivity  of  the  fluid  is  proportional  to  pressure. 
This  situation  has  been  explored  in  Chapter  III. 


APPENDIX  C 
UNCERTAINTY  IN  THE  MEASUREMENT  OF  K. 


The  uncertainty  analysis  presented  here  deals  with 
the  uncertainty  in  the  physical  measurement  of  the  apparent 
thermal  conductivity  coefficient,  K  .   The  apparent  thermal 
conductivity  coefficient,  as  given  by  Equation  42,  is 
a  function  of 

Ka  =  Ka(E,I,6,ro,ri,L,TH,Tc)  (54) 

Equation  54  does  not  explicitly  contain  the  variable  pres- 
sure, P,  under  which  the  insulating  material  was  tested. 
Hence,  the  uncertainty  in  the  measurement  of  K  will  be 
given  with  K  as  the  dependent  variable  and  pressure,  P, 
as  the  independent  variable. 

Before  proceeding  to  calculate  uncertainty  in  the 
measurement  of  K  ,  an  explanation  of  the  difficulty  of  in- 
eluding  pressure  as  an  explicit  variable  in  the  analysis 
is  in  order.   If  K   in  Equation  41  is  replaced  with  the 
sum  of  Equations  19  and  37 

Ka  '  T&T)    C^HT^  d/p3+  W1 

C  n  tt  K   D  N. 

+  s, c  A—  (551 

2  Nh  In  (Dc/y*)  L"J 

then  variables  are  introduced  whose  values  are  not  readily 
available,  such  as  k  ,  d",  n,  N.,  N,  ,  D  .   Using  this 
approach  would  require  assuming  values  for  these  parameters 
and  this  would  negate  the  purpose  of  calculating  uncertainty. 
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The  standard  deviation,  AK  ,  in  the  measurement  of 
K„  is  calculated  with  [37] 


j   3K   2 

*Ka  =  C?  ^    »\02  W 

1     l  z  . 

where 

z  =  independent  variable 

u  •  =  instrument  deviation  in  measuring  z. 

Equation  55  is  based  on  the  assumption  that  the  standard 
deviation  of  the  input  parameter,  z,  is  Gaussian  in  nature. 
Applying  Equation  56  to  Equation  42, 

EI  In  (rQ/ri) 
Ka  =  2ttL  (Th-Tc)  (42) 


and  dividing  the  result  by  K  yields 

3. 


AK  „      -u 

a     Uy,  ?    ux  2       r 

K 
a      -       _  , 

o    "01 


O     ^  0   1V 

u  u 

r-  T„   2     -uT  2  u 

trili(r/r,))  +  2(T-T  )   +  (-T1)  I"         C573 


o'  xi'       XH  *L 

The  maximum  deviation  of  each  of  the  instruments 

used  in  the  measurement  of  the  parameters  in  Equation  54 

and  for  the  pressure  measurement  is  given  below. 
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Parameter 

Temperature 
Current 

Potential 

Length 
Pressure 


Instrument 


Maximum 
Division 


Chromel-Alumel  thermocouple   +  0.75% 

Ammeter-Westinghouse         +  2.5% 
Type  PA- 4 

Weston-Schlumberger  Digital   +  1.0% 
VOM 


Stainless  Steel  Scale 
McLeod  Vacuum  Gauge 


+  1.3xl0-3  ft 
+  3.0% 


Two  sample  problems  to  illustrate  the  range  of  uncer- 
tainty will  be  calculated  for  typical  conditions,  at 
both  high  and  low  temperatures,  encountered  in  the  measure- 
ment of  the  apparent  thermal  conductivity  coefficient 

Sample  1. 

K   =  0.0072  Btu/hrft°F  (Table  5,  line  5) 

TR  =  730°F  (Table  5,  line  5) 

Tc  =  104°F  (Table  5,  line  5) 

E   =  25.0V 

I   =0.55  amperes 

L   =  1.0  ft 

rQ  =  0.163  ft 

ri  =  0.089  ft 


Applying  these  parameters  and  the  instrument  deviations 
from  the  above  list  into  Equation  57  yields 
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AKa   _    rr0.01   x    25.0, 2         r0.025    x   0.55,2         f      0.0013  2 

K  Ll         25.0  J  l  0.55  J  l0.163   x    0.605J 

a 


f      0.0013  .2    +    7    r0-0075    x    750, 2 

l0.089    x   0.605J  l    730    -    104      } 


AK 

-jp  =  0.041 
a 

That  is  AK   =  +  4.1%  of  K   @  P  =  1000  mm  Hg  +  31.   Note 
a  a 

that  in  the   numerator  of  the  first  three  terms  in  the 
above  sample  calculation  that  the  instrument  deviation, 
given  as  a  percentage  of  the  scale  reading,  was  multiplied 
by  the  scaled  value  to  obtain  the  instrument  deviation. 

S  amp 1 e  2 

K   =  0.00407  Btu/hrft°F  (Table  5,  line  8) 

TH  =  452°F  (Table  5,  line  8) 

Tc  =  87°F  (Table  5,  line  8) 

E   =  11.3  V 

I   =0.4  amperes 

L,  r  ,  r-  -  same  as  above 
'   o '   1 

Calculating  as  before,  one  obtains 


^  =    0.041 
Ka 
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or,  AK   =  +  4.1%  6  K   @  P  =  1000  ym  Hg  +  3%. 

3.     ~  3. 

From  the  above  calculations  it  is  apparent  that 
deviations  in  the  measurement  of  K   fall  within  an  accept- 

a. 

able  range. 

Of  the  variables  measured  for  calculating  K  ,  the 
variable  temperature  requires  additional  comments  concern- 
ing methodology  in  the  use  of  the  temperature  transducers 
or  thermocouples.   All  thermocouples  junctions  were  welded 
together  to  form  a  small  bead.   The  thermocouples  used 
to  measure  T„  were  brazed  to  a  copper  surface  and  all 
excess  brazing  material  was  removed.   The  asbestos-covered 
lead  wires  from  the  junction  were  placed  adjacent  to  the 
copper  surface  to  minimize  conduction  losses  from  the 
thermocouple  junction.   The  thermocouple  used  to  measure 
Tf  was  welded  together  and  the  resulting  beaded  junction 
was  ground  away  until  a  thin,  flat  junction  was  obtained. 
This  junction  was  covered  with  a  small  piece  of  ceramic 
insulation  and  then  held  against  a  surface  for  temperature 
measurement.   Efforts  were  made  to  maintain  constant  pressure 
on  the  temperature  probe  when  a  surface  temperature  was 
taken.   This  procedure  yields  more  consistent  data  [38]. 

The  cold  junction  of  the  thermocouple  was  an  ice 
bath.   Even  though  ordinary  tap  water  was  used  in  the  bath, 
a  high  quality  mercury- in-glass  thermometer  showed 
negligible  deviation  from  32°F. 
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